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> ' ABSTRACT 

A comprehensive analysis of 355 high-quality WSRT H I 21 -cm line maps of nearby galaxies 
shows that the properties and incident rate of Damped Lyman a absorption systems (DLAs) 
observed in the spectra of high redshift QSOs are in good agreement with DLAs originating 
in gas disks of galaxies like those in the z « population. Comparison of low z DLA statis- 
tics with the Hi incidence rate and column density distribution f(Nm) for the local galaxy 
, sample shows no evidence for evolution in the integral "cross section density" < na >= l^ 1 

(Z=mean free path between absorbers) below z w 1.5, implying that there is no need for a 
hidden population of galaxies or H I clouds to contribute significantly to the DLA cross sec- 
tion. Compared with z « 4, our data indicates evolution of a factor of two in the comoving 
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O ' density along a line of sight. We find that dN/dz{z = 0) = 0.045 ± 0.006. The idea that 

\ the local galaxy population can explain the DLAs is further strengthened by comparing the 

properties of DLAs and DLA galaxies with the expectations based on our analysis of local 
galaxies. The distribution of luminosities of DLA host galaxies, and of impact parameters be- 
tween QSOs and the centres of DLA galaxies, are in good agreement with what is expected 
from local galaxies. Approximately 87 per cent of low z DLA galaxies are expected to be 
fainter than and 37 per cent have impact parameters less than 1" at z = 0.5. The analy- 
sis shows that some host galaxies with very low impact parameters and low luminosities are 
expected to be missed in optical follow up surveys. The well-known metallicity-luminosity 
relation in galaxies, in combination with metallicity gradients in galaxy disks, cause the ex- 
pected median metallicity of low redshift DLAs to be low (~l/7 solar), which is also in good 
agreement with observations of low z DLAs. We find that f(Njn) can be fitted satisfacto- 
rily with a gamma distribution, a single power-law is not a good fit at the highest column 
densities A^hi > 10 21 cm~ 2 . The vast majority (w 81 per cent) of the Hi gas in the local 
Universe resides in column densities above the classical DLA limit (Nm > 2 x 10 20 cm -2 ), 
with A/hi ~ 10 21 cm -2 dominating the cosmic H I mass density. 
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1 INTRODUCTION 

How does the cosmological mass density of neutral hydrogen (H I) 
gas evolve over the history of the Universe and what sort of galaxies 
are responsible for this evolution? Two completely unrelated ob- 
servational techniques are used to find answers to these questions 
for the present epoch and at earlier cosmic times. At intermediate 
and high redshifts (out to z ~ 6) deep optical and UV spectra of 



background quasars are scrutinised to find high column density H I 
that causes optical depth tll > 1 blueward of the Lyman limit. 
In detecting these absorbing systems, their distance is not a limit- 
ing factor since the detection depends only on the brightness of the 
background source against which the absorber is found. However, 
the nature of the absorbing system is difficult to determine because 
the absorption spectrum only gives information along a very narrow 
sight-line through the system. Despite the large effort that has been 
dedicated to identify and characterise the high redshift absorbers 
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through deep optical imagi ng, very little about their properties is 
yet known (e.g., M0ller et al. 2002, and references therein). 

At z = 0, 21 -cm emission line surveys are used to study 
the high column density Hi. Blind surveys with the Arecibo and 
Parkes radio telescopes have resulted in acc urate me asurements 
of the cosmic Hi mass density Qm(z = )j|Zwaanrtalj |l997t 
iRosenberg & Schneider! 120021 : IZwaan et alJl2003l 120051) . In con- 
trast to the high redshift studies, the 21-cm emission line surveys 
readily result in measurements of the total Hi masses of the de- 
tected objects. Furthermore, after follow-up with an aperture syn- 
thesis instrument, the 3-dimensional H I data cube can be obtained, 
from which a detailed velocity field and column density map can be 
derived. Unfortunately, the Hi 21-cm hyperfine line is very weak, 
which limits the studies to the very local Universe (z < 0.2). The 
column densities that routine 21-cm emission line observations can 
sense, are the same as those of the highest column density Lya 
absorbers known as Damped Lya systems, or DLAs, with column 
densities Ahi > 2 x 10 20 cm -2 . These are the systems that contain 
most of the cosmic H I mass (e.g.. | Wolfe et al. ll986tlLanzetta et alj 
ll99ltlStorrie-Lombardi & Wolfel200Cl) . although there are indica- 
tions that systems with slightly lower column densiti es may con- 
tribute approximately 20 per cent to Qui IPeroux et al. 2005). 

The purpose of this paper is to test whether the results from 
the 2 = 21-cm line surveys and the z > UV and opti- 
cal surveys can be reconciled. This question relates directly to the 
problem of the origin of DLA systems. Traditionally, DLAs were 
thought to arise in large gaseous disk s in the process of evolving 
to present-day spiral galaxies I Wol fe et al. ll9M [T995). This idea 
was supported by the fact that the velocity profiles of unsaturated, 
low-ion metal lines are consistent with rapidly rotating, large thick 
disks (Prochaska &Wolfell99lll998t). Dissenting views do exist, 
perhap s most notably presented by H aehnelt. Steinmetz. & Rauchl 
( 1998), who argued that the velocity structure can alternatively be 
explained by protogalactic clumps coalescing on dark matter haloes 



(see also Ledoux et al. 1998; Khersonskv & Turnshek 1996). In 
this view, the DLA systems are aggregates of dense clouds with 
complex kinematics rather than ordered rotating gas disks. 

Imaging surveys for DLA host galaxies ha ve so far not 
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tifications show that a mixed population of galaxies is responsible 
for the DLA cross section. Semi-analytical and numerical models 
of galaxy formation point to sub-i* galaxies as the major contrib- 
utors to Hi cross section above the DLA limit (e.g., |Kauffniann| 
| l996tlGardner et all200lUNagamine. Springel. & Hern quist 2004; 
lokoshi & Nagashimj2005l) . 

In t his paper we concentr ate primarily on cross-section ar- 
guments. Burbi dge et alJ ll977h started the discussion on whether 
the disks of normal galaxies can explain the incidence rate of ab- 
sorbers. The method consists of multiplying the space density of 
galaxies (measured through the optical luminosity function) with 
the averag e area of the Hi disk above the DLA column den- 
sity limit. iBurbidge et alj ll977h concluded that the mean free 
path for interception with a galactic disk was too large to ac - 
count for the number of observed absorbers. IWolfe et al. I ll986l). 
lLanzetta etai] <199ll). IFvnbo. Moller. & Warrenl <1999h . Ischavd 
(2001), and Chen & Lanzetta 1 2003) use variations of this same 
simple analytical approach for estimates at high and low z. Al- 
though these approaches are useful to gain insight in the problem, 



a much more direct method is to use real 21-cm column density 
maps of local galaxies to evaluate the DLA incidence rate. This 
method was first used bv lRao&BrIggslll993h who used Gaussian 
fits to the radial H I distribution of 27 galaxies as measured with 
the Arecibo Telescope. Much higher resolution H I maps were used 
by IZwaan. Briggs. & Verheiier] l2Q02f) who used WSRT maps of 
a volume-limite d sample of 49 galaxies in the Ursa Major clus- 
ter l Verheiien & Sancisi 2001). A similar s tudy w as conducted 
bv iRvan- Weber. Webster. & Stav elev-Smith 1 2003) based on 35 
galaxies selected from HIPASS iMever et al]|2004l). which were 
observed with the ATCA. Finally, Rosenberg & Schn eider! <2003l) 
used low-resolution (~ 45") VLA maps of their sample of 50 H I- 
selected galaxies. 

With the completion of WHISP (The Westerbork H I survey of 
Spira l and Irregular galaxies. |van der Hulst. van Albada. & SancisH 
2001) we have available a much larger sample of galaxies, with a 
high dynamic range in galaxy properties and observed with a spa- 
tial resolution of ~ 12" x 12"/ sin<5. With this sample, the cross 
section analysis can be done more accurately and in more detail. 
This paper is organised as follows. In the following section we de- 
scribe the details of the WHISP sample, in section 3 we summarise 
surveys for low redshift DLA host galaxies and compile a sample of 
these systems that we use for our comparison with z — 0. We cal- 
culate the H I column density distribution function in section 4 and 
compare the results with the statistics at higher redshifts in section 
5. A comparison of the redshift number density dN/dz between 
2 = and higher redshift is given in section 6. In section 7 we cal- 
culate the probability distribution functions of various parameters 
of cross-section selected galaxy samples and show that these agree 
well with results of low-z DLA host galaxy searches. In Section 8 
we show that the observed metallicities in DLAs are in good agree- 
ment with the expected values if low-2 DLAs arise in the gas disks 
of galaxies typical of the z = population. Section 9 summarises 
the results. We use Ho = 75 kms" 1 Mpc" 1 throughout this paper 
to calculate distance dependent quantities. 



2 THE WHISP SAMPLE 

The incidence rate of high column density Hi in the local 
Universe can best be calculated from 21-cm maps of a large 
sample of galaxies with uniform selection criteria and high 
spatial resolution. Such a sample is available through WHISP 
(The W esterbor k Hi survey of Spiral a nd Irregular galaxies, 
Ivan der Hulst. van Albada. & Sancisil200lh . an H I survey of some 
390 galaxies carried out with the Westerbork Synthesis Radio Tele- 
scope in the periods 1992 to 1998 and 2000 to 2002. The aim of 
this survey is to obtain maps of the distribution and velocity struc- 
ture of the H I in a large sample of galaxies, covering all Hubble 
types from SO to Im and a considerable range in luminosity. This 
survey increases the number of well-analysed H I observations of 
galaxies by an order of magnitude, with the principal goal of inves- 
tigating the systematics of rotation curves and the properties of the 
dark matter haloes they t race. First results on the late type galax- 
ies have been published ISwaters et al]|20 02|) and a study of th e 
early type galaxies is on its way iNoordermeer et al]|2004ll2005h . 
Examples of H I maps of nine galaxies used in our analysis are 
shown in Fig. Q where the top three galaxies have H I masses of 
Mm ~ 1O 8 M , the middle three have M H i ~ 1O 9 M , and the 
bottom three have Mm ~ 10 10 M Q . 

The sample from which WHISP candidates have been se- 
lected consists of galaxies in the Uppsala General Catalogue (UGC) 
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Figure 1. Hi column density maps of six of the WHISP galaxies. The top three galaxies have total Hi masses of 10 8 Mq, the middle three have Mm ~ 
10 9 Mq, and the bottom three have Mm ~ 1O 1O M0. Contour levels are at column densities of 0.2, 0.5, 1, 2, 4, 8 and 16 X 10 21 cm -2 , where the lowest 
contour corresponds to the DLA limit. The beam sizes are indicated in the lower left corners of each panel. The ellipses indicate the orientation parameters of 
the optical images of these galaxies and the extent corresponds to the diameter measured at the 25th B-magnitude isophote. The scale bar in the lower right 
corner corresponds to 5 kpc at the distance of the galaxies. 



of Galaxies iNilson Il973l) that fit the following criteria: 1) B- 
band major axis diameter more extended than 1.2'; 2) declination 
(B1950) north of +20°; and 3) the 21-cm flux density averaged 
over the H I profile exceeds 100 mJy. For observations after the year 
2000 the flux density criterion was relaxed to also include galaxies 
with lower average flux densities. This enabled observation of a 
reasonably large number of galaxies of Hubble type SO - Sb, which 
typically have lower H I fluxes than galaxies of Hubble type Sbc 
and later. For the analysis presented in this paper we made use of 



the data of 355 galaxies out of the total of 391 in the WHISP sam- 
ple, which had been fully reduced at the time this analysis started. 
The distribution over Hubble type and optical luminosity for the 
355 galaxies in our sample is given in table 1. 

A total of about 4500 hours of observing time (more than one 
year of continuous data taking) was used over a period of 10 years 
to observe 391 galaxies. Of these, 281 were observed with the old 
WSRT system between 1992 and 1998. The remaining 1 10 galaxies 
were observed with the upgraded, more sensitive system between 
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Table 1. Properties of the WHISP sample 





L B < L* B /20 


L* B /20 <L B < L* B /5 


L%/5 <L B <L* B 


L B > L* B 


Total 


E/SO 


2 


2 


8 


9 


21 


Sa/Sb 


4 


1 


31 


41 


83 


Sbc/Sc 


7 


14 


37 


52 


110 


Scd/Sd 


19 


15 


17 


4 


55 


Sm/Irr 


43 


31 


7 


5 


86 


Total 


75 


69 


100 


111 


355 



2000 and 2002. The main difference between the two periods is the 
upgrade of the WSRT with cooled receivers on all telescopes and 
a greatly improved bandwidth and correlator capacity. In practice 
this means that the sensitivity of the observations made in 2000 
and thereafter are a factor of ~ 4 better than before (e.g. a typical 
r.m.s. of 0.8 mjy/beam as compared to 3 mjy/beam). The veloc- 
ity resolution for most galaxies is 5 km s~\ The highly inclined 
galaxies in the sample with inclinations typically above 60° have 
been observed at a lower velocity resolution of 20 km s _1 , while 
more face-on galaxies and galaxies with narrow single dish profiles 
(Av ^ 60kms~ ) have been observed with a velocity resolution 
of 2.5 kms" 1 . 

The WSRT observations of the WHISP galaxies have been 
processed following a well described standard procedure consist- 
ing of a (u,v) data processing phase (data editing, calibration and 
Fourier transformation) using the native WSRT reduction pack- 
age NewStar and an image processing phase (continuum sub- 
traction, CLEANing, determination of Hi distributions and ve- 
locity fields using GIPSY. The resulting data cubes and images 
have been archived and a summary is displayed on the web at 
http://www.astro.rug. nl/^whisp A precise description of the re- 
duction procedures can also be found there. 

For each observed galaxy the WHISP pipeline reduction 
recipe produces data cubes with an angular resolution of ~ 12" x 
12"/ sin (5, 30" x 30" and 60" x 60". The maps with different 
resolutions are used to obtain some insight in the effect of spatial 
resolution on the measurement of H I cross sectional areas (see sec- 
tion The median 3-a column density sensitivity limit of the 
highest resolution H I maps is 2.3 x 10 20 cm -2 . At the lower reso- 
lutions (30" and 60") the column density sensitivities improve with 
factors of approximately 3.5 and 10, respectively. To illustrate the 
resolution of our data, we calculate the total number of beams over 
the optical area within the 25th magnitude isophote, and the num- 
ber of beams over the area with an H I column density exceeding 
the DLA limit. These statistics are presented in Fig. [2] The median 
number of beams over the optical area is ~ 40, whereas the median 
number of beams over the H I area is « 220. We note that the to- 
tal number of independent H I column density measurements above 
the DLA limit in our sample is « 140, 000. 

A possible concern about the selection of the WHISP sam- 
ple is that the sample is based on the UGC, which implies that the 
parent sample from which the WHISP galaxies are chosen is op- 
tically selected. Since DLA systems are purely Hi-selected, our 
z = comparison sample of galaxies should ideally also be H I- 
selected to guarantee that all possible optical selection biases are 
ruled out. However, Zwaan 1 2000) has shown that the properties of 
Hl-selected galaxies are not different from those of optically se- 
lected galaxies if compared at the same luminosity. Specifically, he 
showed that the gas richness, optical surface brightness, and scale 
length of optically selected and H I-selected galaxies are indistin- 
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Figure 2. Left: Distribution of number of beams over the optical area within 
D25 ■ Right: Distribution of number of beams over the area with H I column 
density exceeding the DLA limit of 2 x 10 20 cm -2 . 



guishable if subsamples with similar median absolute magnitude 
are compared. Therefore, if the weighting of the sample is prop- 
erly taken into account (using the luminosity function or H I mass 
function) no strong biases should be introduced by using optically 
selected galaxies. As will be discussed in section [4] we use mor- 
phological type-specific H I mass functions to calculate weights for 
individual galaxies, which will ensure that the morphological type 
distribution of our weighted sample at a given H I mass is similar 
to that of an H I-selected galaxy sample. Small differences in the 
optical properties of optically and H I-selected samples as a func- 
tionof Hi mass might remain, but given the findings of IZwaar] 
(2000) and the fact that we apply type-specific weighting, these 
differences are unlikely to bias are results significantly. This is also 
illustrated by the fact that the value for the redshift number den- 
sity dN/dz that we derive in se ction |6| i s consistent with what 
iRvan- Weber. Webster. & Stavelev-SmithTfe003l.l2005al) find based 
on a smaller sample of H I-selected galaxies. 

The field of view of the 21 -cm WSRT observations is defined 
by the diameter of the individual dishes and is approximately 30' 
across (at half power beam width). The WHISP data base contains 
many instances of multiple H I detections within one such field. 
We separated these multiple detections into individual galaxies if 
both Hi detections could be identified with optically identifiable 
galaxies in the NASA/IPAC Extragalactic Database (NED). In most 
cases the multiple detections are known pairs or groups of galax- 
ies. If a second H I detection in the field targeted on a UGC galaxy 
could not be associated with a catalogued galaxy, we did not sep- 
arate this H I detection but treated it as a companion H I cloud to 
the known optical galaxy and therefore included it in the total cross 
section calculation, keeping information about the impact parame- 
ter at which the JVhi occurs. 

All 21 -cm maps of the resulting sample of galaxies were trans- 
formed to column density maps using the well-known equation 
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where Tb is the observed brightness temperature. This equation 
is valid for optically thin emission. For cold, high column density 
(> 10 21 cm~ 2 ) gas clouds, the optically thin approximation might 
break down, which would result in slightly underestimated H I col- 
umn densities. This effect is probably only important for the very 
highest H I column densities, in the range where there are only very 
few DLA measurements. Therefore, it is unlikely to affect our com- 
parison between DLAs and local galaxies. Furthermore, the mea- 
surement of total cross-sectional area above the DLA column den- 
sity limit is not affected at all by H I optical depth effects. 



3 SEARCHES FOR LOW REDSHIFT DLA GALAXIES 

The aim of this paper is to make a detailed comparison between 
the local galaxy population and local DLA galaxies. The difficulty 
in this comparison is that the number of identified DLA galaxies is 
very low, at any redshift. There are now app roximately 600 DLA 
systems known in the Universe IProchaska et all2 005f). only ~ 50 
of which are at z < 1.7. This number is especially small at low 
redshifts, because the Ly-a line is only observable from the ground 
when it is redshifted to z > 1.7. Furthermore, due to the expansion 
of the Universe, the expected number of absorbers along a line of 
sight decrease with decreasing redshift, which naturally leads to a 
scarcity of DLA systems at low z. 

At high z there is a substantial number of systems available, 
but attempts to directly im age their host galaxies have gen erally 
been unsuccessful (see e.g. , Colbert & Malkan 2002; M0ll er et alJ 
l2002t|Prochaska et alfeool) . A few positive identificat ions do ex- 
ist, mostly the result of HST imaging fpy nbo. M0ller, & Warren 
ll999tlWarren et all200ll : lKulkarni et alfeOOCtllVMler et all2002h . 
Possibly the best data to date are available for three objects imaged 
with STIS by M ailer et alj|2002l) . They found emission from three 
DLA galaxies with spectroscopic confirmation and concluded that 
the objects are consistent with being drawn from the population of 
Lyman-break galaxies. 

Although the absolute number of DLAs at low z is small, the 
success rate for finding their host galaxies is better for obvious rea- 
sons: the host galaxies are expected to be brighter and the sepa- 
ration on the sky between the bright QSO and the DLA galaxy is 
likely larger. The first syste matic survey for low z DLA host galax- 
ies was performed bv lLe Brun et alJ |1997L who obtained HST R 
and B-band imaging of seven QSO fields with known DLAs. For 
most of their targets they found likely galaxies, the properties of 
which span a large range from LSB dwarfs to large spirals. Simi- 
lar s earches for single syst ems were done bv lBurbidge et alJ ll996t) 
and lPetitiean et all 11996T) . The DLA at z = 0.656 in the sight- 
line toward 3C33 6 was s tudied extensively bv lSteideletai]ll997l) 
and Bo uche et alJ 1200 lh . but despite the deep HST imaging in the 
optical and Ha, no galaxy has been identified. All these studies 
lacked spectroscopic follow-up, so associations were purely based 
on the proximity of identified galaxies to the QSO sight line, which 
is not unique in many cases. Optical and infrared imaging com- 
bined with spectroscopic follow-up observations was first done by 
iTurnshek et alJ i200lh . who successfully identified two galaxies at 
z = .091 and z — 0.221. Using a similar techniquejLac^ytl 
(2003) also identified two low z DLA galaxies. IChen & La nzetta 
(2003) demonstrated that it is possible to use photometric red- 
shifts to determine the association of optically identified galaxies 



with DLAs, and applied this technique to six DLA systems. Fi- 
nally, Rao et al. (2003) studied four DLA systems and were able 
to confirm two DLA galaxies based on spectroscopic confirmation, 
one using a photometric redshift, and one association was based on 
proximity to the QSO. 

Two attempts have been made to measure directly the cold gas 
contents of low z D LA systems by means of deep 21 -cm emission 
line observations I Kanek aTet alJ 1200 it lLane. Briggs. & S mette 
2000). The problem with these observations is that the 21 -cm line is 
extremely weak so that with present technology surveys are limited 
to the very local (z < 0.2) Universe. Both observations resulted in 
non-detections, which allowed the authors to put upper limits to the 
Hi masses of approximately Mm = 2.3 x 10 9 Mp>, which is one 
third of the H I mass of an L» galaxy {Zwa an et alj2003ll2005h . 

In addition to these surveys for low z DLA galaxies, 
there are two cases of high column density absorption systems 
found in local galaxies. [M iller. K nezek. & Bregmanl 1 19991) and 
iBowen. Tripp. & Jenkinsl 12001) study the absorption line systems 
seen in NGC 4203 and SBS 1543+593, at redshifts of z = 0.004 
and z — 0.009, respectively. For these systems, the galaxy-QSO 
alignment was identified before the absorption was found and in 
one case (NGC 4203) Ly-a absorption has not been seen to date, 
but the structure and strength of observed metal lines resemble 
those of DLAs. 

All together, including the z w systems NGC 4203 and SBS 
1543+593, there are now 20 DLA galaxies known at z < 1. Ta- 
ble [2] summarises the properties of these galaxies. We have been 
fairly generous in assembling this table because not all listed galax- 
ies have been spectroscopically confirmed. Also, some of the sys- 
tems fall just below the classical DLA column density limit of 
log iVm ~ 20.3, but we decided to include these because our 21- 
cm emission line data is also sensitive below this limit. However, 
if in the remainder of this paper when we refer to "DLA column 
densities," we restrict our analysis to log iVm > 20.3 and dis- 
regard the identified systems with lower column densities. Some 
systems in our compilation have been studied by several authors, 
in which case we choose the parameters from the most recent ref- 
erence. The measurements of L/L, are mostly based on B-band 
data, but in some cases these data were not available in which cases 
we used K or i?-band. For the value of L, we adopted the recent 
measu rement of Mg — 51og/i75 = —20 .3 from iNorberg et alJ 
(2002), Mr - 51ogh 75 = - 21.1 from iLin et all 1 1996k and 
M* K - 51og/i 75 = -24.1 from lCole et alJ 1200 ll) . The morpho- 
logical classifications are copied directly from the relevant refer- 
ences and are very heterogeneous in their degree of detail. There- 
fore, we note that the L/L, values and the listed types should be 
treated with caution and both parameters have large uncertainties. 
Our local DLA galaxy sample is defined to have redshifts z < 1, 
which is of course a random choice. However, cutting the sample at 
lower redshift would result in too small a sample to make a mean- 
ingful comparison with our nearby galaxy statistics. The median 
redshift of the DLA galaxy sample is < z >= 0.5. In the analysis 
we ignore any evolutionary effect in the redshift range z = 1 to 
z = 0. This assumption is justified by our finding in section|6| that 
the cross section times comoving space density of DLA systems is 
not evolving in this redshift range. However, an important caveat 
is that this lack of evolution in cross section puts no constraints 
on how gas is distributed in individual systems, but only on the 
total gas content averaged over the whole galaxy population. The 
mean properties of gas disks in galaxies might evolve since z = 1. 
Given the fact that the cosmic star formation rate density drops by 
a factor 8 between z = 1 and z = (e.g.. Hopk insl2004l) . it is pos- 
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sible that evolution of the neutral gas disks actually takes place as 
well. To study the H I mass evolution of galaxies deep observations 
with future instrument such as the Square Kilometer Array (SKA, 
ICarilli & Rawlins J2004) are required. At present no observational 
constraints exist. Apart from the evolution of the H I disks, the op- 
tical properties of DLA host galaxies could evolve betwe en z = 1 
and 2 = as well. However, IChen & Lanzetta 1 2003) find that 
DLA galaxies do not show luminosity evolution between z = 1 
and z — 0, although the constraints are limited by small number 
statistics. 

Finally, one might be concerned about the process of choosing 
which galaxy to attribute to a certain absorber (see discussion in 
IChen & Lanzett all2003l) . In some cases more than one candidate 
is available, in which case the authors are forced to choose which 
one is the most likely absorber. The set of arguments on which 
this choice is based, varies between different studies. We note that 
this introduces some randomness in our compilation presented in 
Tabled 

The conclusion from studying Table|2|is that the sample of low 
z DLA galaxies spans a wide range in galaxy properties, ranging 
from inconspicuous LSB dwarfs to giant spirals and even early type 
galaxies. Obviously, it is not just the luminous, high surface bright- 
ness spiral galaxies that contribute to the H I cross section above the 
DLA threshold, although it is these galaxies that contribute most to 
the total comoving density SIhi- 



4 CALCULATING THE COLUMN DENSITY 
DISTRIBUTION FUNCTION - f(N m ) 

In QSO absorption line studies, the column density distribution 
f{Nm) function is defined such that f(Nm)dNmdX is the num- 
ber of absorbers with column density between TVhi and Nm + 
dNm over an absorption distance interval dX. The analysis in 
this paper derives f(Nm) from the statistics of Hi emission line 
distributions in the local galaxy population in order to predict the 
effectiveness of present day galaxies at accounting for the absorp- 
tion line statistics. The parameter dX assures that f(Nm) counts 
the number of absorbers per co-moving unit of length, but since 
dX/dz — 1 at z = 0, we can replace dX with dz in our calcula- 
tion of /(JVhi). The local f(Nm) can now be calculated from 

J{ m) Ho N m InlOAlogATm ' 

Here, <&(xj) is the space density of objects with property Xj equal 
to that of galaxy i, In reality, the parameter x could be the H I mass 
or the optical luminosity of the galaxies in the sample, so that $(x) 
is the H I mass function or the optical luminosity function of galax- 
ies in the local Universe. Ai (log JVhi) is the area function that de- 
scribes for galaxy i the area in Mpc 2 corresponding to a column 
density in the range log iVm to log Nm + A log Nm- In practice, 
this is simply calculated by summing for each galaxy the number of 
pixels in a certain log Nm range multiplied by the physical area of 
a pixel. The function to(xj) is a weighting function that takes into 
account the varying number of galaxies across the full stretch of x, 
and is calculated by taking the reciprocal of the number of galax- 
ies in the range logx — A/2 to logx + A/2, where A is taken 
to be 0.3. The results are very insensitive to the exact value of A. 
The summation sign denotes a summation over all galaxies in the 
sample. Finally, c/Hq converts the number of systems per Mpc to 
that per unit redshift. Note that dependencies on Ho disappear in 
the final evaluation of f(Nm)- 



The WHISP sample that is used in this analysis is some- 
what biased toward early type galaxies (SO to Sb) because 
these were specifically targeted for the projects described in 
Ivan der Hulst. van Albada. & Sancisl 1200 lh . If we were to apply a 
simple weighting scheme using a luminosity or H I mass function, 
early type galaxies would be given too much weight in the calcula- 
tion of the f(Nm)- To avoid this problem, we choose to app l y type - 
specific Hi mass functions as published by IZwaan et alJ |2003). 
Morphological types are available in LEDA for all but 18 galax- 
ies in the WHISP sample. We divide the WHISP sample into five 
subsets of galaxies: E-S0, Sa-Sb, Sbc-Sc, Scd-Sd, and Sm-Irr. Vi- 
sual inspection of the unclassified galaxies learned that these were 
mostly dwarf irregular or small compact galaxies and hence we put 
these in the Sm/Irr subset. We note that the difference between the 
f(Nm) calculated using the weights based on one Hi mass func- 
tion for all WHISP galaxies and that based on type-specific weight- 
ing is very small. 

Fig. [3] shows our measured column density distribu- 
tion function f(Nm)- We use a binning of logiVm = 
0.1 dex. The error bars indicate la uncertainties and in- 
clude uncertainties in the Hi mass function normalisation from 
Zwaan et al. 1 2003) as well as counting statistics of the WHISP 
sample. The open symbols indicate column densities below 
log Nm = 19.8, which we adopt as our sensitivity limit 
(see Sectio n 14. 1 | for a discussion on this limit). Fo llowing 
| Pei & Falll ll995l). IStorrie-Lombardi. Irwin. & McMahod Jl996h . 
Istorrie-Lombardi & Wolfd feOOOh . and lPeroux et alJ l2003h . we fit 
f(Nm) with a gamma distribution, analogous to the Schechter 
function, which is used for fitting luminosity functions and H I mass 
functions: 

f(Nm) = if* /N^)(Nm/N^)-Pe- N ™? N ™. (3) 

This gives an excellent fit to our data with parameters: j3 — 1.24, 
logA^i = 21.2 cm" 2 , and /* = 0.0193 if we use all data 
points above logAfn > 19.8. The parametric f(Nm) is shown 
as a solid line in Fig. [3] If we restrict the fit to column densities 
above the DLA limit, we find (3 = 1.52, log JVhi = 21.3 cm" 2 , 
and /* = 0.0137, which is shown by the dotted line. We note 
that there is no physical motivation for fitting the /(JVhi) with 
a gamma distribution other than that it provides a reasonable fit. 
Neither is there a physical reaso n to fit f(Nm ) with a power law, 
which is traditionally done (e.g., Tvtler 1987). However, the good 
gamma distribution fit demonstrates that at high H I column den- 
sities (log iVm > 21), the /(JVhi) deviates strongly from a tradi- 
tional single power- law fit. This is con sistent with the form of the 
cut-off identified by Peroux et al. (2003) based on a large sample 
of DLAs and sub-DLAs. 

The total Hi mass density can be calculated by integrating 
over the column density distribution function. When using the 
parametrised gamma distribution, it can be easily seen that the H I 
mass density can be calculated as 

PH i = r(2-/3)/*iV HI mnHo/c, 

where nm is the mass of the hydrogen atom and F is the Euler 
gamma-function. Using the best-fit parameters of our gamma distri- 
bution fit to log iVm > 19.8, we find phi = 6.8 x 10 7 M Q Mpc" 3 , 
in excellent agreement with Z waan et al] 120031) . This may not 
come as a surprise, since we used the H I mass functions from that 
paper to calibrate /(JVhi). Nonetheless, the result of this calcula- 
tion is a good consistency check and assures that the normalisation 
of our f{Nm) is correct (see also lRao & Briggsll993h . 
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Table 2. Properties of z < 1 (sub-)DLA host galaxies. 



QSO 


^abs 


log 7Y H i 
(cm" 2 ) 


b 

(kpc) 


L/L* 


Morphology 


Reference 


Redshift 1 


Ton 1480 


0.004 


20.34 


7.9 


0.44 


early type 


a 


spectro-z 


HS 1543+5921 


0.009 


20.34 


0.4 


0.02 


early type 


b 


spectro-z 


Q0738+313 


0.091 


21.18 


<3.1 


0.08 


LSB 


c 


no 


PKS 0439-433 


0.101 


20.00 


6.8 


1.00 


disk 


d,g,k 


spectro-z 


Q0738+313 


0.221 


20.90 


17.3 


0.17 


dwarf spiral 


eg 


spectro-z 


PKS 0952+179 


0.239 


21.32 


<3.9 


0.02 


dwarf LSB 


e 


photo-z 


PKS 1127-145 


0.313 


21.71 


<5.6 


0.16 


Patchy irr LSB 2 


e,g,l 


spectro-z 


PKS 1229-021 


0.395 


20.75 


6.6 


0.17 


Irr LSB 


f 


no 


Q0809+483 


0.437 


20.80 


8.3 


2.80 


giant Sbc 


f,j,k 


spectro-z 


AO 0235+164 


0.524 


21.70 


12.5 


1.90 


late-type spiral 2 


g,m 


spectro-z 


B2 0827+243 


0.525 


20.30 


29.5 


1.04 


disturbed spiral 


e.g 


spectro-z 


PKS 1629+120 


0.531 


20.69 


14.7 


0.78 


spiral 


e.g 


spectro-z 


LBQS 0058+0155 


0.613 


20.08 


7.5 


0.60 


spiral 


h.k 


spectro-z 


Q 1209+ 107 


0.630 


20.20 


9.7 


1.90 


spiral 


f 


no 


HE 1122-1649 


0.681 


20.45 


23.5 


0.58 


compact 


g 


photo-z 


Q1328+307 


0.692 


21.19 


5.7 


0.67 


LSB 


f 


no 


FBQS 1137+3907C 


0.720 


21.10 


10.3 


0.20 


spiral? 


i 


spectro-z 


FBQS 0051+0041A 


0.740 


20.40 


22.4 


0.80 


spiral? 


i 


spectro-z 


MC 1331+170 


0.744 


21.17 


25.1 


3.00 


edge-on spiral 


f 


no 


PKS 0454+039 


0.860 


20.69 


5.5 


0.60 


compact 


f,j 


no 



References: (a) Miller, Knezek, & Bregmar] jl999l) : fb^ lBowen. Tri pp. & Jenkins 1 2 003): (c^lTurnshek et alj|200ll) : 
(d^ IPetitiean et all Il996h: fel iRao et all <2003t): ffl iLe Bnin et al] ll997t): (e) IChen & Lanzettal 120031): Of) 
Pettini et all feOOfJ) : fi llLacv et all <2003l): (j ) IColbert & Malkanl l2002t) : Ort lChen. Kennicutt. & Rauchl l2005t) : (1) 
Lane et alJll998h : fnfl lBurbidge et alJll99fil) 

1 Indication of whether the redshift of the DLA galaxy has been confirmed by either spectroscopy or a photometric 
redshift measurement 

2 Probably arising in galaxy group 



There are a few interesting features in Fig. [3] First, our mea- 
sured points begin to deviate from the power law end of the gamma 
distribution fit at column densities log Am < 19.8. The first ob- 
vious explanation for this is that in this range of Ami we start 
to loose sensitivity in the 21 -cm maps. There might also be a 
more physical origin of this effect, which is the expected ionisa- 
tion of the outer Hi disks of galaxie s by the m etagalacti c UV- 
background be l ow log Am ~ 19-5 ICorbelli & Salpeten Il993t 
lMalone\<[T993l) . ICorbe lli & Bandied <2002l) argue that if f(Nm) 
is corrected to include H II, it would follow a power-law down to 
Am = 10 17 cm" 2 . Our 21-cm maps lack the sensitivity to study 
this effect in detail. At the high Am end in Fig.[3] the /(Am ) drops 
off exponentially, which is faster than Nm~ 3 above log Am = 
21.6. A fall-off of Am -3 is theoretically expected for randomly 
orie nted gas disk, irrespective of their radial Am distribution 
(see | Milgromlll988tlFan & Peilll993tlZwaan. Verheiien. & Briggsl 
1999). Three independent effects can cause deviation of the Am 
expectation: 1) beam smearing in the 21-cm maps can smooth out 
the very highest column densities; 2) galaxies are not infinitely thin 
(this is assumed in the /(Am) calculation of randomly oriented 
disks); and 3) the Hi gas is not optically thin at the highest col- 
umn densities. In the remainder of this paper, our possible biases 
at the extreme edges of /(A/hi) are not important because we only 
compare our data with the absorption line data in the intermediate 
Am range. Furthermore, the cosmological H I mass density is also 
dominated by the intermediate Am values, as will become clear in 
section[5] 



4.1 The effect of spatial resolution on /(A/hi) 

One concern when comparing column densities measured from 21- 
cm maps with those measured from absorption line systems is the 
large difference in spatial resolution. The median distance of the 
WHISP sample is 20 Mpc, implying that for the highest resolution 
maps the synthesised beam corresponds to « 1.3 kpc. This is al- 
most 10 5 times greater than the diam eter of the optical em ission 
region of a QSO « 0.03 pc, e.g.[Wvithe, A gol. & Flukel2002h . the 
area probed by DLA column density measurements. 

Despite this disparity in resolution, 21-cm emission and Ly- 
man a absorption me asurements generally seem to be in agree- 
ment. iDjckex&iockman 1 1990) noted a consistency between Ly- 
man q absorption measurements towards high latitude stars and 
21-cm emission measurement in the Galaxy. There is only one 
example where a column density from both 21-cm emission and 
DLA absorption has been measured of the same source. The 
DLA galaxy SBS 1543+593 has Lyman a absorption measured at 

2.2 x 10 20 cm" 2 <Bowen. Tripp. & Jenkinj200ll) . whereas the 21- 
cm emission at the position of the background QSO is measured to 
be 5 x 10 20 cm -2 Ichengalu r & Kane kaj2002l) . 

An alternative method to test the effect of resolu- 
tion is to spatially smooth high resolution 21-cm maps 
and measure the resulting column densities and /(A/hi). 
iRvan- Weber. Stavelev-Smith. & WebsteJ l2005bl) used a high res- 
olution 21-cm data cube of the Large Magellanic Cloud and con- 
volved it with circular 2-dimensional Gaussians of various widths 
in the spatial plane. The Am distribution was then re-measured and 
/(A/hi) was re-calculated, based on the new, low resolution column 
densities. Lowering the resolution was found to have a truncating 
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Figure 3. The H I column density distribution function at z = from 21 -cm 
emission line observations. Error bars include uncertainties in the H I mass 
function as well as counting statistics and indicate la uncertainties. Open 
symbols correspond to measurements below our sensitivity limit. The solid 
line is a gamma distribution fit (see to all points above log iVm > 

19.8, the dotted line is the same but restricted to log Am > 20.3. 

effect on the Am distribution, generally decreasing the occurrence 
of the highest column densities. 

For our set of WSRT maps, we can test the effect of lower- 
ing the spatial resolution. This is shown in Fig|4] where we present 
the /(Am) distribution derived from the 30" and 60" maps, along 
with that from the original resolution maps. It can be seen that go- 
ing to lower resolution, the highest column densities are smoothed 
away and this flux appears again as lower column densities. Effec- 
tively, this causes a 'tilt' of /(Am). From this sample it is impos- 
sible to estimate if this tilting would continue if even higher resolu- 
tion maps would be used. However, for our comparison with DLA 
data the effect of resolution is probably unimportant because the 
effects are minimal in the intermediate Am range where we make 
the comparison. 

The redshift number density, dN/dz (see section |6j is 
als o only minimally affected by resolution. The LMC study 
by iRvan- Weber. Stavelev-Smith. & Websted J2005bl) showed that 
dN/dz (at a limit of log Am = 20.3) increases on the order of 
10 per cent when the spatial resolution is changed from 15 pc to 
1.3 kpc, the typical resolution of the WHISP sample. A 10 per 
cent change in the measured dN/dz for the WHISP sample is well 
within the errors quoted in section|6| 

Finally, we comment on our adopted sensitivity limit of 
log Am = 19.8 for the column density distribution function. In 
Section [2] we quoted a formal 3a sensitivity limit of log Am = 
20.3 for the high resolution maps, which implies that below this 
limit we might start to underestimate /(Am). However, Fig. |4| 
shows that /(Am) based on lower resolution data is marginally 
higher in the range 19.8 < log Am < 20.3, and only significantly 
deviates from the high resolution curve for log Am < 19.8. As ex- 
plained above, the higher /(Am) values for lower resolution data 
are not only the result of these data having lower sensitivity lim- 
its and thus picking up lower column density H I, but also because 
small regions of high Am gas are smoothed away and appear again 



Figure 4. The effect of spatial resolution of the 21 -cm line observations on 
the H I column density distribution function at z = 0. Filled circles corre- 
spond to the highest resolution, which is approximately 14", open squares 
and crosses are for the spatially smoothed data, to resolutions of 30" and 
60", respectively. 

as lower column densities. This also explains why small differences 
between /(Am) based on the different resolution data already be- 
come apparent above log Am = 20.3. Therefore, if we were to 
construct a "hybrid" /(Am) using high Am results from high 
resolution data and low Am from low resolution data, we would 
double-count some of the Hi atoms. We choose to simply adopt 
a sensitivity limit of log Am = 19.8, and note that the low end 
(log Am < 20.3) of /(Am) might be slightly underestimated, but 
by not more than 20 per cent. This does not influence any of our 
conclusions presented in this paper. 



5 COMPARISON OF /(Am) WITH HIGH z DATA 

In Fig.|5|we reproduce the /(Am) measurements from our analy- 
sis o f 21-cm maps and co mpare these with data at higher red shifts 
from [Peroux et alJ l2005h. iRao et alj J2005h and IProchaska et alJ 
(2005). The Peroux et al. 1 2005) measurements of /(Am) be- 
low the DLA limit are the result of their new UVES survey 
for "sub-DLAs". We choose here to only plot the points for 
log Am > 19.8, which roughly corresponds to our sensitiv- 
ity cut-off. Their higher column density points stem f r om th e 
combined data o f IStorrie-Lombar di. Irwin. & McMahon| |l99fj). 
IStorrie-Lombardi & Wolfe! <2000h . and IPeroux et alj l200ll). "to- 
gether based on « 100 z > 4 quasars. Thejprochaska et al. 1 2005) 
data are from an automatic search for z > 2.2 DLA systems in 
SDSS-DR3 and also include the results from previous DLA sur- 
veys. We represent their results by both a gamma distribution fit and 
a double power-law fit to the whole data set in the redshift range 
2.2 < z < 5.5. The intermediate redshift points from lRao et alJ 
(2005) are based on Mg II-selected DLA systems. All calculations 
are based on a Sl m = 0.3, S1a = 0.7 cosmology. The surpris- 
ing result from this figure is that there appears to be only very 
mild evolution in the intersection cross section of Hi from redshift 
z ~ 5 to the present. This is a conclusion in stark contrast to ear- 
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Figure 5. Comparison between the H I column density distribution function 
at z = and higher redshifts. The solid points are from our 21-cm emission 
line obse rvations and are the same as in Fig . E] The higher redshift points 
are from IPeroux et all {2005) and R ao et all l2005t) . Horizontal errorbars 
indicate the bin sizes. The curves show the gamma distribution and double 
power-law fits to the SDSS-DR3 results from Prochaska et al. 12005) and 
cover a redshift range of 2.2 < z < 5.5. All results based onaSl m = 0.3 
and S1a = 0.7 cosmology. 



lier works th at claimed strong evoluti on in th e DLA cross section 
(e.g.. lWolfe et al. lll98flLanzetta et aJl99ll) . We will come back 
to this weak evolution in section|6|in which we discuss the redshift 
number density. 

In Fig.|6|we plot the contribution of systems with different 
column densities to the integral Hi mass density pm- At z = 
it is clear that systems with column densities A^m ~ 10 21 cm -2 
dominate the Hi mass density. At higher redshifts (z > 1.7), the 
uncertainties are much larger, but also there it app ears that Am ~ 
10 21 c m" 2 systems contribute most, although the Prochask a et alJ 
1 2005) data suggest that column densities between log Am =20.3 
and 21.3 contribute almost evenly. Any possible differences be- 
tween the results at various redshifts are much more pronounced in 
Fig.|S|than in Fig.[5]because the vertical scale is stretched from nine 
decades to only four decade s. The one point that clearly deviates is 
the highest Am point from iRao et all i2005h at log N m = 21.65. 
This elevated interception rate of high Am Mg II-selected interme- 
diate r edshift DLAs is present in both redshift bins that Ra o et alJ 
120051) distinguish (0.1 < z < 0. 9 and 0.9 < z < 1.7) , and is 
also reported in their earlier work iRao & T urnshek 2000). Fig. [5] 
very clearly demonstrates that this point dominates the f2m mea- 
surement at intermediate redshifts. It is therefore important to un- 
derstand whether the Mg II-based results really indicate that high 
column densities (log Am ~ 21.65) are rare at high redshift, then 
indeed become more abundant at intermediate redshifts (0.1 < z < 
1.7) and subsequently evolve to the scarce numbers at z — 0. Al- 
ternatively, the high /(A/hi) point might be a result of yet uniden- 
tified selection effects introduced by the Mg II-selection, although 
presently there are no indications that support this I Pero ux et alJ 
l2004lRao et alj2005t) . 

To further investigate the relative contribution of different col- 
umn densities to the H I mass density, we plot in Fig. the cu- 



Figure 6. The Hi mass density contained in systems of different col- 
umn density per decade of Am- Solid circles are for z = and are 
from our data, othe r symbols are for higher redshifts and are derived from 
|Perouxetal]|2005t) and lRao et aJ] l2005t). The curves are the converted fits 
to the SDDS-DR3 results from lProchaska et"ai]<2005l) . 

mulative H I mass density as a function of column density. The 
solid curve is based on our measured f(Nm) points, not on the 
gamma distribution fit to the data. From this curve alone, we 
would conclude that at z = the fractional mass in systems with 
log Am > 20.3 (classical DLAs) is 86 per cent. However, the con- 
tribution of gas with column densities below our sensitivity limit of 
approximately log A^hi = 19.8 is very uncertain. Unfortunately, 
only selective regions have been imaged to low column density 
limits, but no large scale 21-cm emission line surveys that reach 
sensitivities below log Am ~ 19.8 are available yet. For exam- 
ple, the M31 environment has been imaged with the GBT and the 
WSR T at resolutions between 50 pc and 1 1 kpc iBraun & ThilkeJ 
12004 and references therein), covering the column density range 
17 < logA^Hi < 22. However, the Hi emission in this region is 
mostly due to discrete High Velocity Clouds (HVCs) physically 
associated to M31. These surveys may not be representative of 
all low column density Hi gas at z = 0, and do not provide 
good constraints on the shape of /(Am). To obtain a rough es- 
timate of the contribution of gas below our sensitivity limit, we 
extrapolate our measured f( Nm) below l og Am = 19.8 with 
f(N m ) oc Am -1 - 5 (per e.g., lTvtleJl987l) and find that the frac- 
tional mass in systems with log Am > 20.3 is ~ 81 per cent, 
and hence only ~ 19 per cent of the Hi atoms in the local Uni- 
verUniverse are in sub-DLA column densities. For illustrative pur- 
poses, we also show as a shaded area the cumulative distribution 
functions based on the gamma distribution fits to f{Nm) as pre- 
sented in Figure [3] Here, the fit to log Am > 19.8 corresponds 
to the upper boundary and the fit to log Am > 20.3 to the lower 
boundary to this area. Obviously, simply extrapolating the gamma 
distribution fits to f(Nm) below the DLA limit results in a severe 
overestimation of the number of H I atoms in sub-DLA systems. 

The higher r edshift curves are based on the data of 
|PerouxetaiT l2005). With their new column density measurements 
of individual sub-DLAs, these autho rs improved the earlier esti- 
mates of f(Nm) presented in lPeroux et alJ i2003l) and now show 
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Figure 7. Cumulative H I mass distribution in systems with column density 
> TVhi ■ The solid line is based on our data alone, the dashed line includes 
a correction for H I gas with column densities below our sensitivity limit 
and indicates that ss 81 per cent of the integral mass in Hi is in classical 
damped systems with log Nm > 20.3 cm -2 . The shaded area indicates 
the distribution based on the gamma distribution fits to our f(Nm) (Fig- 
urel3l. where the lower boundary corresponds to the fit to log A^hi > 20.3 
and the upper boundary to the fit to log N^j > 19.8. The long-dashed and 
dash-dotted lines are from lPeroux et alj 12005) and are for redshift ranges 
1.8 < z < 3.5 and z > 3.5, respectively. 

that the fractional contribution of sub-DLAs to Qui is approxi- 
mately 20 per cent at all redshifts z > 1.8. These curves also indi- 
cate that at high redshifts the highest column densities are relatively 
rarer, as can also be seen in Figure[5] Prochaska et al. 1 2005) esti- 
mate the fractional contribution of sub-DLA systems to the integral 
mass density by combining their meas ured f(Nm) with th e red- 
shift number density of sub-DLAs from lPeroux et alj 1200 it) . They 
conclude that the sub-DLAs (or super-LLSs in their terminology) 
contribute between 20 and 50 per cent over the redshift range 2.2 to 
5.5, but emphasise that their single power-law fit to f(Nm) might 
even underestimate this fraction. Obviously, the importance of sub- 
DLA systems at high redshift remains a topic needing further clar- 
ification in the future. 

The fact that from our high resolution 21 -cm maps we find that 
~ 81 per cent of the H I mass density at z = is in column den- 
sities above the DLA limit, implies that Qui measurements from 
blind 21-cm surveys overestimate £Iola(z = 0) by only w 23 per 
cent. When comparing measurements of the atomic gas mass den- 
sity at high and low redshifts, it is important that the z — value be 
corrected with a factor 0.81 when compared to measurements from 
DLAs. On the other hand, when at higher redshift the sub-DLA 
contribution is taken into account (e.g., Peroux et al. 2005), no cor- 
rection of the z — value is required. We refer to lProchaska et al] 
(2005) for various definitions of fi eas and their recommended us- 



age. 



6 THE REDSHIFT NUMBER DENSITY AT z = 

The redshift number density dN/dz for column densities larger 
than Nui can be calculated from the z = f(Nm) distribution by 



Figure 8. Integral redshift number density dN/ dz as a function of H I col- 
umn density cut-off. For log Nm > 20.3 cm -2 we find dN/dz = 0.045, 
dropping rapidly for higher column densities. 



summing over all column densities larger than Nui : 



dN/dz 



N'>N H 



c 

Tf 



E 



£.$(xQu)(x0A(logiV£ T ) 
A log AT' 



(4) 



(5) 



where definitions of $, w, and Ai are equal to those in section|4] 
In Fig. [8] we show dN/dz as a function of JVhi- For column den- 
sities in excess of the DLA limit of logTVui = 20.3 we find 
dN/dz = 0.045 ± 0.006, where the la uncertainty is domi- 
nated by the uncertainty in the parameters of the Hi mass func- 
tion that is used to calibrate /(JVhi). For larger column densi- 
ties dN/dz drops rapidly, systems with logA r Hi > 21 contribute 
only 20 per cent to the total DLA cross section. Our measure- 
ment of dN/dz translates into a mean cross section density of 
< no >= (1.13 ± 0.15) x 10" 5 Mpc -1 , or to a mean free path 
between absorbers of I —< no > -1 = 88 ± 12 Gpc. 

The value of dN/dz(z = 0) = 0.045 agrees very well with 
our previous measurement of dN/dz = 0.042 ± 0.015, based on 
a mu c h smaller sample of galaxies jZwaan. Briggs. & Verheiienl 
2002). Rosenberg & Schneider 1 2003) used Hl-selected galax- 
ies to find dN/dz = 0.053 ± 0.013, also in good agreement 
with our value. Their slightly higher value is probably the re- 
sult of the steeper Hi m ass function slope that these authors use 
to cali brate their /"(iVwTl. lRvan- Weber. Webster. & Stavelev-Smithl 
(2003) also used a sample of Hl-selected galaxies observed with 
the ATCA and applied the same HI mass function normalisation 
as in the present analysis, to find dN/dz = 0. 0461° ^ (see also 
iRvan-We be^ Webste^ & St avelev-Smithl2005a|) . An earlier calcu- 
lation bv lRao^^BriggT^993h resulted in a much lower value of 
dN/dz = 0.015. This discrepancy arises partly because these 
authors limited their analysis to large, optically bright galaxies, 
whereas the newer values take into account dwarf and low sur- 
face brightness (LSB) galaxies, and partly because their calcula- 
tion is based on the Gaussian luminosity function from iTammannl 
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1 1985), which has a lower normalisation than the most recent Hi 
mass function measurements. 

It has been noted before that the area of the H I disk of galax- 
ies above the DLA limit correlates very tightly to the H I mass (e.g., 
lBroeilJl992tlRosenberg & Schneideill2003l) . We can make use of 
this correlation to obtain a alternative measurement of dN/dz. 
Suppose that the projected Hi area A correlates with Mm as 
A oc /3Mm, and A* is the projected area of an Mm* galaxy, then 

dN/dz = J e(Mm)A(Mm)dMm (6) 

= -^A*9*r(l + a + f3), (7) 

where 6** and a are the normalisation and the low-mass power-law 
slope of the H I mass function, respectively. For our WHISP sam- 
ple we find that A* — 992 kpc 2 . Using this value, the lZwaan et alJ 
(2003) Hi mass function parameters (a = —1.30 and 8* = 
0.0086 Mpc" 3 ), and (3 = 1, we find that dN/dz = 0.044. Quite 
surprisingly, this simplistic approach yields a result very close to 
our measurement of dN/dz = 0.045 based on a much more de- 
tailed analysis. 

6.1 Evolution of dN/dz 

We now have a robust measurement of the redshift number den- 
sity at z = 0. How does this value compare to DLA measure- 
ments at higher redshifts? In Fig. [5] we show the combined re- 
sults of high and low redshift dN/dz measurements from dif- 
feren t surveys. At z = 0, we plot our poin t s as well as those 
from iRvan- Weber. Webster. & Stavelev-Smithl 120031 12005 al) and 
lRosenber^^^dmeideJl2003l) . All these local values are based on 
analyses of Hi 21-cm maps of galaxies, have small error bars and 
are consist ent with each other. The points at 2 ~ 0.6 and z w 1.2 
are from lRao et alJ 1200 5) based on their z < 1.65 HST survey of 
Mg II selected systems from the SDSS-EDR. Their total sample of 
z < 1.65 DLAs is now 41, which is a considerable improvement 
over th e original sample of 16 DLAs presented in lRao & Turnshekl 
( 2000). The underlying assumption that the incidence of DLA sys- 
tems can be derived from Mg II surveys relies on the empirical fact 
that all DLA systems show Mg II absorption, whereas the reverse 
is not true. The known incidence of Mg II systems can therefore 
be used to bootstra p the DLA statistics. A similar technique was 
used bv lChurchilll ^20011) , who found dN/dz = 0.08t ;^ at a 
median redshift of z = 0.05. This number is based on HST data on 
four Mg II absorbers. The high redshift points from Prochaska et al. 
(2005) are the result of an automatic search for z > 2.2 DLA sys- 
tems in SDSS-DR3 and also include the results from previous DLA 
surveys as summar ised in IStorrie-Lombardi & Wolfel fcOOd) and 
|Perouxetal]l2003h . 

Also shown in Fig. |5] by a dash ed line is the fit dN /dz = 
0.055(1 + z) 1 11 from St orrie-Lombardi & Wolfel fcOOfJ) . which 
agrees reasonably well with the z — value. It is often quoted 
in the literature that this fit indicates "no intrinsic evolution in the 
product of space density and cross-section" of damped absorbers, 
which is another way of saying that the number of systems per co- 
moving unit of length does not evolve. This statement is based on 
the fact that for a q = Universe, dX/dz = 1 + z £s (1 + z) 111 . 
However, for a modern non-zero A Universe, dX/dz is given by 

£ = - {1 + Z)2 (8) 

dz Vn M (l + z) 3 - (!2m + - 1)(1 + zf + n A 

which starts to deviate significantly from the qo — prediction for 
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Figure 9. Number density of DLAs per unit redshift as function of redshift. 
All points at z > 0.1 are from UV and optical surveys for DLA systems or 
Mgll systems (see text). The points at z ~ are from techniques similar 
to that presented in this work. The dashed line shows the best fit to dN/ dz 
from Storrie-Lombardi & Wolfe 1 200Q). The solid line represents "no evo- 
lution in the product of cross section and co-moving space density" for a 
cosmology with f2 m = 0.3 and Ca = 0.7. This line is scaled vertically so 
as to fit our z = point. 

2 > 1. For no evolution in the number of systems per co-moving 
unit of length, we expect 

dN = /dJVN dX 
dz \ dz ) z= o dz 

The solid line in Fig. [5] represents this run of dN/dz as function 
of z for an Q m = 0.3 and Oa = 0.7 Universe. We normalise this 
line to our dN/dz(z = 0) = 0.045 measurement. With respect to 
this fi m = 0.3, Qa = 0.7 prediction, there is only weak evolution 
in the comoving incidence rate from z ~ 4 to the present time. 
Between z ~ 1.5 and 2 = there is no evidence for evolution 
at all, between 2 ~ 4 and 2 = 0, the evolution in the comoving 
incidence rate is approximately a factor 2. 

This conclusion contrasts with previous claims that the local 
galaxy population cannot explain the DLA incidence rate. For ex- 
ample, estimates of t he cross sec tion to DLA absorpt ion in local 
galaxy disks bv lWolfe et al.Ul986t) and lLanzetta et aljjl99ll) indi- 
cated that there should be evolution of at least a fact or of 2 to 4 (de- 
pendin g on the value of go)- Likewise, the analysis o flRao &*B riggs 
( 1993) pointed toward strong evolution in the incidence rate since 
2 = 2.5. Part of the reason for our conclusion being different from 
older works, is the change in cosmological parameters. For no evo- 
lution in the comoving number density in a go = 0.5 cosmology, 
the allowed change in dN/ dz is much smaller than that for a mod- 
ern fl m = 0.3, Q,a = 0.7 cosmology. The other reason is that the 
galaxy luminosity functions that were used for older dN/dz cal- 
culations had a lower normalisation than the more recent estimates 
from large scale optical and 21-cm surveys. 

We emphasise that the normalisation of f(Nm) and hence 
of dN/dz is completely independent of the WHISP galaxy sam- 
ple, but instead depends only on the Hi mass functions derived 
from HIPASS. This is a large scale blind 21-cm emission line sur- 
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vey covering the whole southern hemisphere and the redshift range 
z = to z = 0.04. Most of the weight to the normalisation of the 
Hi mass function comes from galaxies around z = 0.01. Optical 
and infrared surveys have shown that a large angular area around 
the southern Galactic pole is underdense by fa 25 per cent e xtend- 
ing out to z = 0.1 jFrith et aljfcoOSt iBusswell et alJl2004l . This 
area occupies approximately one third of the HIPASS sky cover- 
age, but it is not clear whether the underdensity extends to even 
larger regions in the southern hemisphere. In any case, our mea- 
sured f(Nm) and dN/ dz might be underestimated by up to 25 per 
cent due to this local galaxy deficiency. This implies that the evolu- 
tion in the comoving incidence rate is perhaps even slightly weaker 
than portrayed in Fig.|5| 

The lack of evolution in the comoving incidence rate since 
redshift z ~ 1.5 implies that the average Hi cross section above 
the DLA limit has not changed significantly over half the age of the 
Universe. At present, it is difficult to ascertain whether this should 
be interpreted as no evolution in the DLA population - meaning no 
change in the space density and size of DLA absorbing systems - 
or as a combined effect where an evolving space density is com- 
pensated by a change in mean absorber size. 



7 EXPECTED PROPERTIES OF LOW-REDSHIFT DLA 
HOST GALAXIES 

In this section we calculate the probability distribution functions 
of the expected properties of galaxies responsible for high column 
density H I absorption at z = 0. To this end, we define a quantity J\f 
as the 'volume density of cross sectional area' for different column 
density cut-offs. J\f is calculated as 



«=-1.3 



E(y,x)$(x)dx, 



(10) 



where E(y,x) is the cross-sectional area of Hi in Mpc 2 above a 
certain column density cut-off as a function of galaxy properties x 
and y, and $(x) is the space density of galaxies in Mpc - as a 
function of x. Similarly to the calculation of f(Nui), the parame- 
ter x has elements A/hi and L, such that $ is again the H I mass 
function or the optical luminosity function that is used to calculate 
real space densities. The vector y could be any galaxy property, 
such as optical surface brightness or morphological type. For ex- 
ample, if x = Mhi and y = (i, we calculate AT(fJ,), the volume 
density of cross-sectional area as function of surface brightness /i, 
using the H I mass function to calculate space densities. In the re- 
mainder of this paper we will use x = Mm, so that the H I mass 
function is used for calibration. We find that the conclusions would 
not change significantly if we were to use the optical luminosity 
function instead. 

Another way of looking at N is that it defines the number of 
systems above a certain column density limit that would be encoun- 
tered along a random 1 Mpc path through the z = Universe. Put 
differently, J\f can be written as dN/dz x Ho/c, where dN/dz, 
the number of systems per unit redshift, is a familiar quantity in 
QSO absorption line studies. The run of M as a function of galaxy 
property x represents the probability distribution of x of galaxies 
responsible for absorption above a certain column density. 

In the following, we present this probability distribution for 
four different column density limits, logA?Hi > 19.8, > 20.3, 
> 20.8, and > 21.3. The distribution for log N m = 20.3, which 
is the classical DLA limit, is always shown as a thick line. In table 
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Figure 10. The expected distribution of absolute B-band magnitudes of 
high column density H I absorbing systems. The lines plus errorbars show 
the product of cross sectional area and space density, which translates to 
the number of expected absorbers per Mpc per magnitude. The right axis 
shows the corresponding number of absorbers per unit redshift dN/dz. The 
different lines correspond to different column density limits, as indicated 
by the labels. The thick line corresponds to the classical DLA limit of 
log TVjji > 20.3. The smaller panels show the effect of changing the slope 
of the low-mass slope of the H I mass function that is used to calculate the 
normalisation. 

13 we tabulate the relative contributions of galaxies of different lu- 
minosities and H I masses to the total DLA cross section. In table|4] 
the mean and median properties of expected DLA host galaxies are 
summarized. 



7.1 Luminosities 

In Figure fTol we show M as a function of B-band absolute mag- 
nitude. As described in the preceding section, this distribution is 
calculated by multiplying the Hi mass function with the cross- 
sectional area above a certain column density cut-off, and bin- 
ning in absolute magnitude. The two lower panels in Figure [Tol 
show the effect of changing the slope of the HI mass function 
that is used to calibrate the distribution. The luminosity measure- 
ments are taken from the RC3, or when these are not available, 
from LEDA, which gives B-band magnitudes transforme d to the 
RC3 s ystem. Since we use the L» measurement from lNorberg et all 
which is in the bj system, we convert our magnitudes to bj 
using bj — Brc3 + 0.185 as derived bv lLiske et al. I2003J), and 
hence Figure fTol is approximately in the bj system. What is imme- 
diately obvious from this plot is that the probability distribution is 
not strongly peaked around L» galaxies. Rather, for column den- 
sities above the DLA limit, the distribution is almost flat between 
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Mb w — 15 and Mb ~ —20. The consequence of this is that if an 
H I column density Nm > 10 20 ' 3 cm -2 is encountered somewhere 
in the local Universe, the probability that this gas is associated with 
an L„/50 galaxy is only slightly lower than for association with 
an L* galaxy. More specifically, 87 per cent of the DLA cross sec- 
tion is in sub-L* galaxies and 45 per cent of the cross section is in 
galaxies with L < L./10. These numbers agree very well with the 
luminosity distribution of z < 1 DLA host galaxies. Taking into 
account the three non-detections of DLA host galaxies and assum- 
ing that these are <C L, , we find that 80 per cent of the z < 1 DLA 
galaxies is sub-L». The median absolute magnitude of a z = 
DLA galaxy is expected to be M B = -18.1 (~ L«/7), with 68 
per cent in the range —15.5 < Mb < —20.0, whereas the mean 
luminosity is L»/2. 5. 

By studying a sample of low-z DLA galaxies, iRao et all 
(2003) also conclude that low luminosity galaxie s dominate the 
Hi cross section, whereas IChen & Lanzetta (2003) claim that lu- 
minous galaxies can explain most of the DLA systems and that a 
contribution by dwarfs is not necessary. The origin of this apparent 
disagreement lies in the de finition of a dwarf galaxy. Using the def- 
inition of Rao et all j2003) that all sub-L» are dwarfs, these galax- 
ies would indeed dominate the DLA cross section. However, using 
the more stringent definition of dwarf galaxies being fainter than 
L*/10, these systems only contribute approximately 45 per cent of 
the cross-section. 

What can also be seen from Fig. llQI is that the highest column 
densities (logiVHi > 21.3) are largely associated with the most 
luminous galaxies: the probability distribution for the highest col- 
umn densities is more peaked around i* . Small cross sections of 
very high column density gas can still account for high H I masses, 
which explains why the total H I mass density is dominated by L» 
galaxies. 

7.2 H I masses 

In Figure 1111 we show JV per decade of H I mass as a function 
of log H I mass. This figure shows largely the same behaviour as 
Figure flol but the distribution is somewhat more peaked around 
Mj^j galaxies (log M^/Mq — 9.8). The contribution from sub- 
Mj^j galaxies to the total DLA cross section is still 81 per cent, 
that of galaxies with H I masses lower tha n 10 9 M(T) is 31 per cent. 
These figu res agree very well with those of Rcwmber g &_SchneideJ 
120031) and lRvan- Weber. Webster. & Stavelev-Smithl <2003l) . based 
on samples of galaxies observed with lower spatial resolution. 
The median Hi mass of a z = DLA galaxy is expected to 
be logA/m/M© = 9.3, with 68 per cent in the range 8.5 < 
log A/hi/Mq < 9.8. Again we see that the highest column den- 
sities are preferentially associated with the most massive galaxies. 
We will come back to this fact in section l7~5l 



7.3 Surface brightness and Hubble type 

The probability distribution function of ^ifj 5 , the mean B-band sur- 
face brightness within th e 25th mag arcsec -2 , is given in F ig. 1 121 
Simila r to the finding of iRvan- W eber. Web ster. & Stavelev- Smith 
( 2003), we find that the cross section is dominated by galaxies with 
p?B m the range 23 to 24 mag arcsec -2 . For reference, for our sam- 
ple we find the median value of pf| for i» galaxies is 23.6 mag 
arcsec -2 . At the bright end the distribution drops rapidly, showing 
that galaxies with high surface brightnesses contain a small frac- 
tion of the cross section. Toward dimmer galaxies, the distribution 
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Figure 11. The expected distribution of the H I masses of high column 
density H I absorbing systems. The lines plus errorbars show the product 
of cross sectional area and space density, which translates to the num- 
ber of expected absorbers per Mpc per decade of H I mass. The right axis 
shows the corresponding number of absorbers per unit redshift dN/dz.The 
different lines correspond to different column density limits, as indicated 
by the labels. The thick line corresponds to the classical DLA limit of 
log iVjji > 20.3. The smaller panels show the effect of changing the slope 
of the low-mass slope of the H I mass function that is used to calculate the 
normalisation. 

drops off slower. We find that 53 per cent of the DLA cross section 
is in galaxies dimmer than /i^ 5 = 23.6 mag arcsec -2 . However, for 
8 per cent of the galaxies in our sample no measurement of ^fj 5 is 
available. Assuming that these galaxies have no measurement be- 
cause of their low surface brightness, we are biased against LSB 
galaxies. Including these in our lowest p?B bins increases the frac- 
tion of cross section in galaxies dimmer than 23.6 mag arcsec -2 to 
64 per cent. 

The measurement of optical surface brightness that we have 
used tends to understate the contribution of LSB galaxies to the 
DLA cross section. The reason for this is that if the surface bright- 
ness is measured within a fixed isophote, the radius of this isophote 
shrinks if the surface brightness decreases. (j,b is therefore always 
measured over the central brightest part of a galaxy, which naturally 
decreases the dynamic range in surface brightness measurements. 
A better alternative would be /tiff , the effective surface brightness 
(the average surface brightness within the half-light radius), but un- 
fortunately this parameter is only available for 40 per cent of the 
galaxies in the WHISP sample. For illustrative purposes, we can 
assign values of fi C B to those galaxies that have no measurements, 
by using the surface brightness-luminosity relation o bserved for 
nearby galaxies, for example bv lCross & Driven 12002). For galax- 
ies without ^b we simply take the measurement of the galaxy 
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Figure 12. The expected distribution of the B-band surface brightness 
within the 25th mag isophote of host galaxies of high column density H I 
absorbing systems. The lines plus errorbars show the product of cross sec- 
tional area and space density, which translates to the number of expected 
absorbers per Mpc per magnitude. The right axis shows the corresponding 
number of absorbers per unit redshift dN/ dz.The different lines correspond 
to different column density limits, as indicated by the labels. The thick line 
corresponds to the classical DLA limit of log Nm > 20.3. 
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Figure 13. The expected distribution of the morphological types of host 
galaxies of high column density H I absorbing systems. The lines plus er- 
rorbars show the product of cross sectional area and space density, which 
translates to the number of expected absorbers per Mpc. The right axis 
shows the corresponding number of absorbers per unit redshift dN/ ciz.The 
different lines correspond to different column density limits, as indicated 
by the labels. The thick line corresponds to the classical DLA limit of 
log A/m > 20.3. 



closest in absolute magnitude in our sample. Thus, we find that 
the probability distribution of fi c g is much flatter at the LSB end 
than that of /Mg. Specifically, we find that 71 per cent of the cross 
section is in galaxies dimmer than fi c g = 22.0 mag arcsec -2 , 
which corresponds to the peak of the distribution for L, galaxies. 
If an LSB galaxy is defined as having a surface brightness fainter 
than 1.5 mag arcsec - below this value, we find that 44 per cent 
of the cross section is in LS B galaxies, in accord with the find- 
ings of iMinchin et alJ 120041) . The fractional contribution of LSB 
galaxies to the DLA cross section is larger than their contribution 
to the H I mass density because their typical H I mass densities are 
lower than those observed in high surface brightness galaxies (e.g., 
Ide Blok & McGaugblll997l). This is su pported by the galaxy for- 
mation models of Mo. Ma o. & White! Il998t) . which indicate that 
cross-section-selected samples are weighted towards galaxies with 
high angular momentum (i.e., LSB galaxies). In conclusion, our 
data is not ideal to make firm statements about the contribution of 
LSB galaxies to the DLA cross-section. However, using the infor- 
mation we have available we can state that galaxies with surface 
brightness dimmer than that of a typical L„ galaxy make up at least 
half of the cross section. 

Turning now to the Hubble types, we see in Fig. ll3l tb.at late- 
type galaxies are preponderant in the distribution of DLA cross sec- 
tion. Earlier types are not negligible and the probability of identi- 
fying an SO galaxy with a DLA is approximately a third of finding 
an Sc galaxy. 



Table 3. Relative contribution to DLA cross section from dif- 
ferent galaxies 



Quantity 


< L./10 


< L./5 


< L« 


< 2L» 


L 


0.45 


0.58 


0.87 


0.96 


A/hi 


0.22 


0.37 


0.81 


0.96 



contributed by each element dNmdb on a fine grid in the Nm-b 
plane. We again used the type-specific Hi mass functions to as- 
sign weights to each individual galaxy (see section[4j- To increase 
the signal-to-noise in the figure, we smoothed the probability dis- 
tribution with a Gaussian filter, which results in a final resolution 
of a — 0.13 dex in the iVm direction and a = 1.3 kpc in the b 
direction. The contour levels are chosen at 10, 30, 50, 70, and 90 
per cent of the maximum value. The shaded area in the figure indi- 
cates the column density region that corresponds to DLA column 
densities. 

A few interesting features can be readily seen in Figure Il4l 
The lowest contour illustrates that the highest Hi column den- 
sities are very rarely seen at large galactocentric radii: there ap- 
pears to be a strong correlation between Nm and the maximum 
radius at which this is observed. The observational fact that the 

Table 4. Expected properties of low-z DLAs 



7.4 Impact parameters 

Figure fl4l shows the probability distribution of H I cross section in 
the Nm-b plane, where b is the impact parameter in kpc from the 
position of the observed column density to the centre of the galaxy. 
The WHISP sample was used to calculate the cross sectional area 



Quantity 


median 


mean 


logarithmic mean 


M B 

log A/hi (Mq) 
b (kpc) 


-18.ll?;S 

9 3+ - 5 

'•°-s.o 


-19.2 
9.5 
10.6 


-17.7 
9.0 
7.0 
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Figure 14. The two-dimensional probability distribution of H I cross section 
in the H I column density - impact parameter plane. The contours are calcu- 
lated from our H I 21 -cm maps and are drawn at 10, 30, 50, 70 and 90 per 
cent of maximum. The points are from DLA galaxy searches in the z < 1 
Universe. The shaded area shows the column density region corresponding 
to classical damped systems (log TVjii > 20.3) 

Hi distribution in galaxies often shows a central depression can 
be seen by the compression of contours near 6 = kpc. Finally, 
this plot shows that the largest concentration of H I cross section in 
galaxies in the local Universe is in column densities in the range 
20.3 < log TVhi < 20.8 and impact parameters b < 7 kpc. 

The points in Figure^] are the pairs of iVm and b measure- 
ments from the literature sample of low z DLA galaxies from Ta- 
ble|2llflow2DLAs are drawn from the same population of galax- 
ies as those in our sample from the local Universe, we would expect 
the points to show the same distribution as the contours. Unfor- 
tunately, the statistics are too poor to calculate contours from the 
DLA galaxy data, but there seems to a qualitative agreement be- 
tween the two data sets. 

A more straightforward comparison between the DLA galax- 
ies and local galaxy properties is presented in Figure Il5l which 
shows the conditional probability distribution of impact parameter 
b as a function of column density iVm- This figure is calculated 
by normalising at each value of JVhi in Fi sure Il4l the distribution 
function of cross section as function of b to the peak of the distri- 
bution. The solid line shows the median b value as a function of 
TYm, the other lines show the 10th, 25th, 75th, 90th, and 99th per- 
centiles, respectively. The literature values are again overplotted as 
points. Both in the z — data and in the DLA da ta, there is a weak 
correlation between JVhi and fo. lRao et alJi20 03') also noted the ex- 
istence of this relation. The agreement between the points and the 
contours is remarkably good: 50 per cent of the literature values 
are within the 25th and 75th percentiles, and 75 per cent are within 
the 10th and 90th percentile contours. Based on our analysis, the 
expected median impact parameter of log iVm > 20.3 systems is 
7.6 kpc, whereas the median impact parameter of identified z < 1 
DLA galaxies is 8.3 kpc. Although we are limited by small num- 
ber statistics, a comparison between the contours and the points 
suggests that the observed number of very low b systems is lower 
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Figure 15. Conditional probability of impact parameter as a function of H I 
column density. The thick solid line shows the median impact parameter as 
calculated from our 21-cm emission line maps. The other lines show the 
10th, 25th, 75th, 90th, and 99th percentiles, as indicated by the labels. The 
points are from DLA galaxy searches in the z < 1 Universe. 

than expected. These identifications are the ones most likely to be 
missed due to the proximity of a the bright background QSO. Al- 
ternatively, if dust obscuration in DLA host galaxies is important, 
these low b sight lines are expected to be the first to drop out of a 
flux limited quasar sample. 

A more detailed view of the distribution of impact parameters 
is presented in Figure fl6l which shows the predicted probability 
distribution of b for different column density cut-offs. The peak 
of the distribution for DLA column densities is at b ~ 5 kpc and 
the distribution drops rapidly toward higher values of b. The inset 
shows the probability distribution of b above the DLA limit on a lin- 
ear scale. For higher column densities, the probability distributions 
drop off even more rapidly, which again shows that it is extremely 
unlikely to encounter a high JVhi at a large separation from the cen- 
tre of a galaxy. For DLA column densities, we find that 60 per cent 
of the host galaxies are expected at impact parameters < 10 kpc 
and 32 per cent at < 5 kpc. Assuming no evolution in the proper- 
ties of galaxy's gas disk, these numbers imply that 37 per cent of 
the impact parameters are expected to be less than 1" for systems 
at z = 0.5 and 48 per cent less than 1" at z = 1. These numbers 
illustrate that very high spatial resolution imaging programs are re- 
quired to successfully identify a typical DLA galaxy at z ~ 1. 

The intrinsic assumption in this comparison is that the low z 
DLA galaxy sample is a fair cross-section selected sample. In re- 
ality the sample is a compilation of many surveys, using different 
selection techniques, and different resolutions and wave bands to 
image the galaxy. Keeping this limitation in mind, we conclude 
that the measured impact parameters and column densities of low z 
DLA galaxies is in agreement with the hypothesis that DLA galax- 
ies can be explained by the local galaxy population. 

For completeness, we show in Figure lTTl the expected proba- 
bility distribution function of impact parameter b at various column 
densities. Note the difference with the previous diagrams, where 
we plotted probability distribution functions above certain column 
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Figure 16. The probability distribution of impact parameter b between the 
background QSO and the centre of a galaxy giving rise to Hi absorp- 
tion. The lines show the product of cross sectional area and space density, 
which translates to the number of expected absorbers per Mpc per kpc. The 
right axis shows the corresponding number of absorbers per unit redshift 
dN/dz.The different lines correspond to different column density limits, 
as indicated by the labels. The thick line corresponds to the classical DLA 
limit of log Nfti > 20.3. The inset shows the probability distribution of b 
for column densities log TVhi > 20.3 on a linear vertical scale. 



1 1 I I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 

\70% 




og N HI [cm- 2 ] 



Figure 18. Probability distribution of optical luminosities of galaxies giving 
rise to H I absorption with column density Nm and impact parameter b. 
The lines indicate the probability that the host galaxy is more luminous 
than an L* galaxy, and represent 5, 10, 30, 50 (thick line) and 70 per cent 
probabilities, from bottom to top. The points are from DLA galaxy searches 
and are the same as in Fig. ll4l and Fig. 1151 but here the symbol size represent 
intrinsic luminosity of the host galaxy. Red points indicate galaxies brighter 
than L„, blue points are fainter than L». 
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Figure 17. Left: Normalised probability distribution of impact parameter b 
for different column densities. For H I column densities near the DLA limit 
of 2 X 10 20 cm -2 , the most likely impact parameter to the host galaxy is 
5 hyg 1 kpc. Right: Normalised probability distribution of H I column den- 
sity for different impact parameters b. 

density cut-offs. The right panel is the probability distribution func- 
tion of column density at various values of b. 



7.5 Linking DLA parameters to galaxy properties 

We conclude our comparisons between DLA systems and local 
galaxies by looking at the combined probability distribution func- 
tion of column density, impact parameter and host galaxy luminos- 
ity. In Fig. ll8l we again show the Nm and b measurements from 
low-z DLA galaxies from the literature, but this time the symbol 
size reflects the luminosity of the galaxies such that the symbol 
area scales in direct proportion to L/L t . The contours represent 
the probabilities that the combined measurement of Nni and b is 
expected for a galaxy with luminosity L > L» . For example, the 
thick solid line divides the diagram in two regions, above this line 



most host galaxies would be more luminous than L„, below this 
line most would be less luminous than L, . Apparently, the most lu- 
minous galaxies are most likely associated with high column den- 
sity DLAs, at large impact parameters from the background QSOs. 

The probability contours can be directly compared to the dis- 
tribution of symbol sizes. Although the scatter is large, we see a 
general agreement between the points and the contours: below the 
5 per cent line mostly low luminosity galaxies are found, and on 
the other hand, the two galaxies above the 70 per cent line are both 
galaxies with high luminosities. A similar diagram is presented in 
Fig. 1191 but here contours indicate the probability that the host 
galaxy has an Hi mass in excess of M m . The conclusions from 
this figure are very similar as those from Figure fT8l 



7.6 Comparison to models and simulations of DLAs 

In order to compare our findings with numbers calculated in models 
of galaxy formation, we first transform the B-band luminosities of 
our WHISP galaxies into rotational velocities via the Tully-Fisher 
relation. We choose here t o use the Tully-Fisher relation deter- 
mined by Mever et al. (2005), based on the HIPASS sample. Fig- 
ure |20| shows the cumulative distribution of DLA redshift number 
density for galaxies with different rotational velocities Kirc We 
present the cumulative distribution because this representation is 
normally used in publications based on cosmological simulations, 
and can therefore be directly compared to th ose. 

On the basis of semi-analytical models Okoshi & Nag ashimal 



1 In principle, a measurement of rotational velocity is available directly 
from the WHISP data, but at present this analysis has not been completed 
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Figure 19. Probability distribution of H I masses of galaxies giving rise to 
H I absorption with column density TVjji and impact parameter b. The lines 
indicate the probability that the host galaxy has an H I mass in excess of 
Mm*, and represent 5, 10, 30, 50 (thick line) and 70 per cent probabilities, 
from bottom to top. The points are the same as in Fig. 1181 
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Figure 20. The cumulative redshift number density of DLAs in galaxies 
of different rotational velocities. The different lines correspond to different 
column density limits, the thick line is for log Nm > 20.3. 

(2005) find that the average virial velocity of a 2 = DLA 
is Vvir ~ 90 kms -1 and the average luminosity is 0.05L,. In 
their models, the fraction of DLA hosts in galaxies fainter than 
L*/10 is 98 per cent, whereas we find that this fraction is only 

for the full sample. Since the observed scatter in the Tully-Fisher relation is 
very small, the conclusions will not change by using this approximation. 



41 per cent. Furthermore, they find that the typical impact param- 
eter is 3 kpc, much smaller that our median value of 7.8 kpc. To 
make their resul ts agree with ob servations of low-z DLA galax- 
ies, lOkoshi & Nagashimal ^2005) propose that the masking effect 
where the DLA galaxies are contaminated by the point spread func- 
tion of the QSO, hinders the identification of 60 per cent to 90 
per cent of DLA galaxies with small impact parameters. The SPH 
simulations of Naeamine, Sprin geL & Hernquistl 120041) also sug- 
gest that the relative contribution of low mass galaxies to the DLA 
cross section is very large: their cumulative contribution of red- 
shift number density as a function of Kirc is steeper than what 
we find in Figure 1201 These authors note that the z — re- 
sults should be taken with caution because the mass resolutions 
that these results are based on are low and higher resolution SPH 
simulations are probably required to arrive at mor e precise re - 
sults at z = 0. Oth e r work, such a s that of |Gardner_etalJ j200l|), 
iMo. Mao. & White! < 19981 Il999h . lHaehnelt. Steinmetz. & Rauchl 
( 2000), and Mailer et al. ( 2001) mostly concentrate on the high red- 
shift (z ~ 3) DLA population and cannot be compared directly to 
our work, but also point to sub-L, galaxies as the major contribu- 
tors to the DLA cross section. In conclusion, semi-analytical mod- 
els and cosmological simulations generally over-predict the redshift 
number density contribution o f low mass systems. A notable ex - 
ception is the simple models of lBoissier. Peroux. & Pettiml 120031) . 
which show that the peak and the median of the cross section dis- 
tribution lies around L» galaxies. These models under-predict the 
importance of low mass systems, both in comparison with our re- 
sults and in comparison with low z DLA host galaxies. 

From Figure|20|we find that the mean rotational velocity V c i TC 
of a z = DLA is 111 kms" 1 (the lo g-weighted mean V c irc is 
95 kms -1 ). Note that the rotational velocities we measure are the 
peak velocities of the galaxies' rotation curves, whereas in simu- 
lations galaxies are normally characterised by their virial velocity 
V^,i r . The ratio between these two values depends on the concentra- 
tion index o f the dark halo, but a typical value is V c irc ~ 1.4Kdr 
(see B ullock et all200ll) . which implies that the log-weighted mean 
virial velocity of a z = DLA is approximately 70 kms -1 . Us- 
ing the relatio n betw een V c irc and the virial mass M V [ T given by 
iBullock et alJ 1200 ll) . we find that the mean total mass of a z — 
DLA would be « 1.5 x 10 n M Q . 



8 METAL ABUNDANCES OF LOW REDSHIFT DLAS 

In the previous section we have presented evidence that DLAs arise 
in the gas disks of galaxies such as those in the z — population. 
Measuring metallicities in DLAs therefore should provide infor- 
mation on abundances of the interstellar matter in galaxies. Deter- 
mining DLA metallicities as a function of redshift will probe the 
history of metal production in galaxies over cosmic time. How- 
ever, the metallicities typically measured in DLAs are low (around 
1/13 solar), which is often taken as an indication that DLAs do 
not trace the general galaxy population, for which a mass-weighted 
mean metallicity of near-solar is expected at low redshift (see e.g., 
IPettini et alll997^ . In this section we present a more detailed anal- 
ysis of the expected metallicities of DLAs, under the assumption 
that they arise in the gas disks of normal present day galaxies. We 
have no direct measurements of metal abundances for the complete 
sample of WHISP galaxies. However, we can make use of metal- 
licity studies in other galaxies to statistically assign metallicities to 
our sample. Although this approach might introduce some degree 
of uncertainty in the results, it will help in understanding the ob- 



18 M. A. Zwaan et al. 



served metal abundance measurements in DLA systems. It should 
be kept in mind that the abundances of local galaxies are measured 
from emission lines arising in the photo-ionised gas, while the mea- 
surements in DLAs are from absorpt ion lines in the neutral gas. 
However, recent work bv lSchulte-Ladbeck et alJ j 2005 ) shows that 
at least in one well-studied nearby galaxy the emission and absorp- 
tion measurements of the same a-elements are consistent. 

8.1 The expected metallicities of DLAs at z = 

We take two different approaches to assigning metallicities to 
our sample of galaxies. The first approach is to adopt the well- 
established metallicity-luminosity (Z — L) relation for local galax- 
ies and apply that to ou r sample. We choose to adopt the relation 
given bv lCiarnettll2002l) for oxygen abundances: 

log (O/H) = -0.16A/ S - 6.4. (11) 

The results do not change significantly if instead we use the more 
recent Z — L relation derived from SDSS imaging and spectroscopy 
of 53,000 galaxies iTremonti et all2004 . To convert the log(0/H) 
values to solar values we ado pt the solar abundance of 12 + 
log (O/H) = 8.66 from lAsplund et alj<2004) . The effect of abun- 
dance gradients as a function of g alactocentric distance is taken 
into ac count by using the result from lFerguson. Gallagher. & Wvsd 
( 1998), who found a mean gradient of [O/H] of —0.09 dexkpc -1 
along the major axes of spiral galaxies. We assume that this gra- 
dient is equal for all galaxies and that it does not vary as a func- 
tion of radius. We apply a simple geometrical correction to correct 
for inclined disks. Furthermore, we assume that Eq. I1 II applies to 
the mean abundance measurement within Rib, which agrees with 
Garnet 3 <2002h . For each galaxy in our sample we scale the offset 
of the gradient such that the mean abundance measurement within 
R25 fits the metallicity-luminosity relation. This approach intrinsi- 
cally assumes that the [O/H] abundance gradients are independent 
of galaxy luminosity. To test the effect of this assumption on the 
calculation of expected metallicities of DLAs, we also take a sec- 
ond approach in which we adopt an [O/H] gradient that has both 
the intercept and slope varying with galaxy absolute magnitude. We 
fit stra ight lines to the relations found bv lVila-Costas & Edmunds 
1 1992) (their Figure 1 1) to find the varying slope and intercept. The 
two approaches depend on completely independent observations of 
metal abundances in local galaxies. 

Every pixel in our 21 -cm maps is assigned [O/H] values, based 
on the two procedures. Applying a similar method to that set out in 
section|7| we can now calculate the expected [O/H] distribution for 
cross-section selected samples, for different column density limits. 
The results are shown in Figure I2TI where the solid lines refer to 
the first approach of fixed gradients, and the dashed lines refer to 
the varying gradients. 

Although the distributions peak at slightly different locations, 
it is obvious that the global trends are not strongly dependent on 
our assumptions on abundance gradients in disks. The main con- 
clusion is that the metallicity distribution for H I column densities 
logTVHi > 20.3 peaks around [0/H]=— 1 to —0.7, much lower 
than the mean value of an L, galaxy of [0/H]~ 0. Furthermore, 
we see a strong correlation between the peak of the [O/H] proba- 
bility distribution and the H I column density limit: for H I column 
densities in excess of logA^i = 21.3, the expected peak of the 
distribution is at [0/H]=— 0.2. This correlation is expected, since 
our assumed radial abundance gradients impose a relation between 
iVm and [O/H]. Taking into account uncertainties in the Hi mass 
function, in the L — Z relation and in the abundance gradients in 
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Figure 21. The probability distribution of oxygen abundance [O/H] of H I 
absorbers. The lines show the product of cross sectional area and space 
density, which translates to the number of expected absorbers per Mpc per 
decade. The right axis shows the corresponding number of absorbers per 
unit redshift dN/dz. The solid lines refer to the approach of assuming fixed 
[O/H] gradients and the the dashed lines refer to varying gradients (see text). 
The four different lines for each approach correspond to different column 
density limits, as indicated by the labels. The thick lines corresponds to the 
classical DLA limit of log N m > 20.3. 

galaxies, we adopt a value of [0/H]= —0.85 ± 0.2 as a representa- 
tive value for the median cross-section-weighted abundance of H I 
gas in the local Universe above the DLA column density limit. 

8.2 Comparison to low-z DLA metallicity measurements 

It is interesting to compare this expected abundance of DLAs at 
z = to higher redshift values to see whether the measure- 
ments can be brought into agreement. The first evidence of evo- 
lution in DLA metallicities over cosmic time was presented by 
Kulkarni & Fall (2002), contrary to previous claims of no evolu- 
tion (e.g., Pettini et al. 1997, 1999). This results was further estab- 
lished by an analysis of 125 DLA metallic ity measurem ents over 
the redshift range 0.5 < z < 5 by Prochaska et al. 1 2003), who re- 
ported significant evolution o f —0.26 ± 0.07 dex p er unit redshift 
in the mean metallicity [M/H], Kulkar ni et al J 120051) compiled new 
metallicity measurements in four DLAs at z < 0.52 and combined 
with literature measurements this study also finds evol ution in 
[M/H] of -0.18 ± 0.07 dex per unit redshift. The lProchaska et alJ 
(2003) metallicity measurements are mostly based on a-element 
abundances, whereas the Kulkarni et al. 1 2005) data is mostly based 
on Zn abundances. We refer to lK^iikarniet^l 12005) for a discus- 
sion on why Zn is an appropriate choice as a metallicity indicator. 

Before we compare our z — results to those for DLAs at 
higher redshifts, we discuss the different measurements of 'mean 
metallicity' that we can use to make the comparison. The first is to 
calculate the global interstellar metallicity Z, defined as 

2_ J Nmf{Nm)Z(Nm)dN m 

J Nmf(Nm)dN m ' ^ 
The reason that this equation gives a true measurement of the 
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interstellar metallicity is that absorption lines are cross-section- 
selected, which implies that a correct weighting of different re- 
gions in abs orbing systems is autom atically taken into account, as 
is stressed by Kulkarni & Fall 1 2002). In practice, in DLA studies Z 
is normally calculated by taking the average over the metallicities 
of individual absorption-line systems weighted by their H I column 
densities. For our low z data we can determine Z by simply taking 
the mass-weighted average of all mean metallicities of the WHISP 
galaxies: 

/ e(M m )M m Z{M m )dM m 
J e(M m )M m dM m ' 

where Q(Mm) is the Hi mass function and Z(Mui) is the 
mean metallicity of a galaxy with Hi mass Mm. We can cal- 
culate Z(Adm) by taking the column density-weighted metallic- 
ity of each WHISP galaxy, but only counting column densities 
above the DLA limit. Applying this method we find values of 
Z — 0.58^0 and Z = 0.35^0 fo r the fi xed gradient and the 
variable gradient me t hod, r espectively. iKulkarni & Fall 120021) and 
iFukugita & Peeblesl j2004) applied similar techniques, but inte- 
grated over the optical luminosity function and found values of 
Z « 0.8Z Q and Z = 0.83 ± O.25Z . The reason for these values 
being higher than ours is that these studies use one global metal- 
licity measurement for each galaxy (which involves some radial 
averaging over the inner parts of galaxies), while our analysis takes 
into account radial abundance gradients over the whole H I gas disk. 
For the mean mass-weighted metallicity of H I gas with log Nm > 
20.3 at z = we adopt the value of log Z/Z @ = -0.35 ± 0.2, 
where the errors again represent all uncertainties in arriving at this 
result. 

The second measurement of metallicity is the cross-section- 
weighted average, which we will refer to as Z. For DLAs, this 
measurement is simply the unweighted mean of the metallicities 
and represents a 'typical' metallicity of DLAs. For our z = 
sample, we can make use of the calculations presented in Fig. 12 II 
which shows the cross-section-weighted [O/H] distribution of high 
column density Hi gas in the local Universe. We find values of 
Z = 0.36^0 and Z = 0.21Z© for the fixed gradient and the 
variable gradient method, respectively, and hence adopt a value of 
log Z/Zq = -0.55 ± 0.2, for the 'typical' metallicity of z = 
DLAs. 

The third method stems directly from the way we presented 
the probability distribution of metallicity in Figure flTI As dis- 
cussed above, we derive from this the cross-section-weighted me- 
dian metallicity, for which we adopted log Z/Zq — —0.85 ± 0.2. 

In Figure l22l we show the metallicities in DLAs as a func- 
tion of redshift , extra cted from the combined compilations by 
IProchaska et alj <2003l) and lKulkarni et alJ <2005l) . In amalgamat- 
ing the two data sets we gave preferen ce to Zn m easurements for 
systems that occurred in both lists. The lKulkarni et all (2005) data 
include many Zn measurements at low z, but a large fraction of 
those are upper limits. In order to take these limits into account ap- 
propriately, we applied the Kaplan-Meier estimator for randomly 
censored data sets, as implemented in the survival analysis package 
ASURV iFeigelson & Nelsonl 19851) . In six different redshift bins 
(approximately the same as those defined by Kulkarni et al. l2005t) 
we calculate the mean A^m-weighted (Z), the mean cross-section- 
weighted (Z) and the median metallicities, as indicated by solid 
circles, triangles and stars, respectively. The errorbars are la sta- 
tistical uncertainties as given by the Kaplan-Meier estimator. The 
long-dashed, short-dashed and dotted lines show the linear least 
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Figure 22. The metallicity [M/H] measurements in DLAs as a function of 
redshift. Data at redshifts z > are taken from IProchaska etafl fc003l) and 
IKulkarni et alj 120051) . Survival analysis has been used to calculate mean 
and median values because the data include many upper limits. The symbols 
refer to column density-weighted mean (filled circles), the cross-section- 
weighted mean (filled triangles), and median values (filled squares) in six 
redshift bins. The errorbars indicate 1 <r statistical uncertainties. The lines 
are least squares fits to the z > data points. The points at z = are from 
our analysis as described in the text. The open symbols at z ~ 0.6 show 
the effect of excluding the X-ray metallicity measurement in Q0235+164. 
The squares and Mangles are offset horizontally by —0.1 dex and 0.1 dex, 
respectively. 

squares fits to each set of measurements. The larger symbols at 
z = show our estimates of metallicities in the local Universe. 

For all measurements of metallicity a slope between —0.25 
and —0.3 dex p er unit redshif t can b e se en, in agreement wit h 
the findings of IProchaska et alJ 120031) and lKulkarni et alj fc005l) . 
on whose compilations of data this figure is based. This increase 
in metallicity is consistent between the highest redshift points at 
z ~ 4 and the lowest points at z ~ 0.6. If we assume that DLAs 
arise in the gas disks of galaxies, where continued star formation 
causes a release of metals into the ISM, it seems reasonable to 
assume that this evolution in metallicity persists down to z = 0. 
Note also that most models of cosmic chemical evolution predict a 
nearly exponential increase in metallicity as a function of redshift 
(c.f.. lKulkarni & Fall 2002). Given this assumption, we find from 
Figure l22l that the expected metallicity of DLAs at z — is in 
excellent agreement with our estimates from the local galaxy pop- 
ulation: the z = extrapolations of all three lines are within the 
uncertaintie s of our estimat es. 

IKulkarni et al] 120051) have argued that the Nui -weighted 
DLA metallicities do not rise up to near-solar values at z = 0. 
These authors find a DLA metal abundance of —0.79 ± O.I8Z0 
at z = u sing survival analysis, a factor of two lower than the 
IProchaska et alj {200j|) extrapolation. This discrepancy lies mostly 
in the exclusion of the x-ray absorption based m etallicity of the 
low-z system AO 0235+164 jjunkkarinen et alJl2004T) . This sys- 
tem causes the large errorbars on the z ~ 0.6 points in Figure [22] 
The open symbols in Fig. 1221 shows the effect of excluding AO 
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0235+164. The mean and median [M/H] are not affected much, 
but the Nm -weighted value drops by more than a factor of two 
because of the high iVm measured in this system. Hence, if we 
only used the UV absorption measurements, we would find that 
our iVm-weighted Z in galaxies is a factor two higher than the ex- 
trapolated value from DLAs. Different authors disagree on whether 
it is fair to include this s ystem (c.f., IChen. Kennicutt. & Rauchl 
l2005t Kulkarni et al. 2005), but at least it illustrates how poor the 
statistics are on low-z abundance measurements. The lowest red- 
shift bin includes 17 measurements (including AO 0235+164), of 
which five are limits. For illustration, if we calculate log Z/Zq 
at z = from 17 random sightlines through our galaxy sample, 
we find that the statistical error on this value is —0.2 + 0.35 dex, 
demonstrating that metallicity measurements based on small sam- 
ples ofr2LAs_are_jmaTOWably_marked by large uncertainties (see 
also lChen. Kennicutt. & Raucrfeooll) . 

Interestingly, we see that for local galaxies the Nm -weighted 
Z is higher than the mean Z. Although the statistics are poor, it 
appears that for z > 1.5 DLAs this is not the case. The origin 
of this difference might lie in the dust obscuration bias in DLAs: 
DLAs with the highest Hi column densities and highest metal- 
licities could be missed in magnitude-lim ited surveys as a conse- 
quence of their own extinction (see e.g. lFaU^^Pej[^89^ HJ93 ; 

[Elliso n et al. l200lLlMurphv & Liskd2004llVladilo & Perouxl2005b 
Wild & Hewettl2005t for a continued discussion on this issue). If, 
like at z = 0, there is an intrinsic positive correlation between iVm 
and metallicity, this bias would have a much higher effect on the 
iVm-weighted Z than on the mean Z of DLAs, causing the two 
measurements to be very similar in magnitude. In our local galaxy 
sample such a bias would not exist. The fact that we see a factor of 
two difference between the two values at z — and not for z > 1.5 
DLAs therefore is suggestive of dust obscuration introducing biases 
in DLA samples. 

Finally, we comment on the scatter in the observed metallicity 
measurements in DLAs. From the probability distribution of [O/H] 
in Fig. 1211 we find that the expected la scatter in cross section- 
selected metallicity measuremen ts is 0.7 dex . A similar value is es- 
timated bv lChen. Kennicutt. & Raucr] f2005). This is slightly larger 
than the typical scatter observed in metallicity measurements at 
high redshift of ~ 0.5 dex, illustrating that metallicity gradients in 
galaxy disks can easily explain the observed scatter in DLA metal- 
licity measurements. 

In conclusion, we find that the cross-section-weighted mean 
metallicity of local galaxies as well as the mass-weighted metal- 
licity are in very good agreement with the hypothesis that DLAs 
arise in the H I disks of galaxies, provided that the metallicities in 
DLAs continue to evolve since z — 0.5 with a rate similar to that 
observed between z — 0.5 and = 4. 

Before ending t his section, we wish to comment briefly on 
the work of IChen. Kennicutt. & Rauchl 120051) . who presented sim- 
ilar conclusions based on Monte-Carlo simulations. These authors 
use average radial H I profiles for three types of galaxies and as- 
sume that the metallicity gradients in these galaxies are similar to 
what is found for their sample of DLAs. Their analysis does not 
take into account that galaxies over a large range in absolute mag- 
nitude contribute to the DLA cross section (see section l7~Ti . and 
effectively, only L„ galaxies are considered. Therefore, the L — Z 
relation is not taken into account in this analysis. Furthermore, in 
their analysis, the mean metallicity is measured over the whole 
gas disk out to the radius where 21-cm emission is normally de- 
tected, and not just over the region log iVm > 20.3. Since there 
is an intrinsic relation between iVm and metallicity, this causes an 



underestimation of the mean Z for DLAs. This explains why the 
IChen. Kennicutt. & Rauchl l2005t) estimate of the Nm -weighted 
mean Z — O.SZq is lower than our estimate. 



9 CONCLUSIONS 

In this paper we tested the hypothesis that DLA absorption lines 
observed in the spectra of background QSOs arise in gas disks of 
galaxies like those in the z — population. Since the H I column 
densities seen in DLA systems (jV H l > 2 x 10 20 ) are the same 
as those routinely observed in 21-cm emission line studies of local 
galaxies, we can make use of these observations to test the hypothe- 
sis. Thus, we used a sample of 355 high quality WSRT 21-cm emis- 
sion line maps to calculate in detail the expected column density 
distribution function, the redshift number density, and the expected 
probability distribution functions of different galaxy parameters of 
the low redshift DLAs. We summarise the conclusions as follows: 

1) The local galaxy population can explain the incidence rate 
of low redshift DLAs. There appears to be no evolution in the 'cross 
section times space density' or the mean free path between ab- 
sorbers from z ~ 1.5 to z = 0. Between the highest redshifts at 
which DLAs are found (z ~ 4 — 5) and the present time, the evo- 
lution in comoving incidence rate is only approximately a factor of 
two. We find that dN/dz(z = 0) = 0.045 ± 0.006. 

2) Based on the local galaxy population it is expected that the 
DLA cross section is dominated by sub-i* galaxies (87 per cent). 
This agrees with the statistics of identified DLA host galaxies at 
low redshifts. 50 per cent of the low redshift DLAs should arise 
in galaxies with H I masses less than M m = 6 x 10 9 M Q . The 
median z = DLA arises in a L*/7 galaxy with an H I mass of 
2 x 1O 9 M . 

3) The distribution of impact parameters and column densities 
agrees very well between local galaxies and low-z DLA galaxies. 
The median impact parameter between the line of sight to a QSO 
and the centre of the galaxy giving rise to a DLA is 7.8 kpc. For 
systems at z — 0.5 (z — 1) we expect that 37 per cent (48 per 
cent) have impact parameters less than 1". These findings support 
indications that optical surveys for DLA host galaxies miss identi- 
fications at very small impact parameters, because of the brightness 
of the QSO or because of blending due to too low spatial resolution. 
If obscuration of background QSOs by dust in DLA galaxies is im- 
portant, this might also have the strongest effect at small impact 
parameters. 

4) We combine our data set with the well-established 
luminosity-metallicity relation of galaxies and observed metallicity 
gradients in galaxy disks to estimate the expected metallicity dis- 
tribution of low- z DLAs. We find that the expected median metal- 
licity of z — DLAs is approximately 1/7 solar, in good agree- 
ment with observations of metal lines in DLAs. The mean mass- 
weighted metallicity of the interstellar matter in local galaxies 
above the DLA limit is approximately half solar. This is consistent 
with extrapolations from higher redshift measurements, although 
the z — extrapolated value has large uncertainties given the poor 
statistics from DLAs with redshifts approximately z ~ 1.5. 

5) The column density distribution function /(JVhi) in the lo- 
cal Universe can be fitted satisfactorily with a gamma distribution. 
A single power law is not a good fit. There is remarkably little evo- 
lution in the shape of f{Nm) from high z to the present. 

6) Most (~ 81 per cent) of the cosmological mass density in 
Hi at z = is locked up in column densities above the classical 
DLA limit of Nui > 2 x 10 20 cm -2 , the rest is mostly in column 
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densities just below this limit. The fraction is consistent over the 
redshift range z ~ 5 to z — 0. 
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